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Preface 

There  are  a  number  of  options  for  modulation  formats  in  optical  fiber  and  free-space  lasercom  systems. 
There  are  a  seemingly  even  greater  number  of  options  for  receiver  architectures.  Although  most  of  these 
possibilities  and  combinations  have  been  analyzed  in  journal  publications  and  textbooks,  there  does  not 
seem  to  be  a  single  place  to  look  for  all  relevant  results  concerning  these  options.  Furthermore, 
comparisons  of  the  options  do  not  seem  to  have  been  methodically  tabulated. 

In  this  work,  then,  we  will  present,  describe,  analyze,  and  compare  the  theoretical  communications 
performance  metrics  of  the  most  common  optical  modulation  formats  paired  with  the  basic  classes  of 
receiver  types. 

Before  we  start,  however,  we  must  warn  the  reader  that  this  material  is  not  presented  in  the  style  of  a 
standard  textbook,  with  the  concepts  painstakingly  developed  and  derived  from  first  principles.  The  style 
we  have  taken  throughout  is  to  present  the  required  mathematical  formulas  and  show  the  performance 
curves,  but  to  give  only  the  sketchiest  of  derivations,  if  any.  Similarly,  detailed  references  are  given  in 
only  some  cases  although  the  more  general  textbook  references  listed  may  help  to  till  in  many  missing 
details.  It  is  hoped  that  the  specific  topics  discussed,  the  formulas  presented,  and  the  several  references 
given,  along  with  some  simple  library  and  web  search  skills,  will  lead  interested  readers  to  as  much  extra 
development  and  detail  as  they  find  interesting  or  useful. 

Despite  what  this  report  is  not ,  though,  we  hope  that,  as  a  compendium,  it  will  find  utility  for  the 
communications  engineer  and  analyst.  We  have  also  tried  to  sprinkle  throughout  comments,  observations, 
and  opinions  that  are  not  generally  found  in  technical  treatises.  We  hope  these  are  useful  to  the  reader. 
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1.  Introduction 

The  most  common  transmission  formats  used  in  optical  communications  systems  include  on-off  keying 
(OOK),  phase  shift  keying  (PSK),  and  the  so-called  orthogonal  modulations — pulse  position  modulation 
(PPM),  frequency  shift  keying  (FSK),  and  orthogonal  polarization  shift  keying  (POLSK).  (Of  course, 
optical  signaling  has  at  least  as  many  other  variants  as  classical  radio  and  electronic  signaling,  but  we  will 
constrain  ourselves  in  this  report  to  the  listed  formats,  which  have  been  used  traditionally.)  For  each  of 
these,  we  will  examine  several  receiver  architectures,  including  preamplified,  heterodyne,  homodyne, 
photon  counting,  and  optimum  quantum.  We  will  examine  both  coherent  and  incoherent  receivers  when 
applicable.  Metrics  to  be  examined  include  the  bit  and  symbol  error  rates  of  uncoded  systems  and  the 
channel  capacity  and  computation  cutoff  rate  of  coded  systems.  For  these,  we  will  examine  both  hard- 
decision  and  maximum-likelihood  (soft-decision)  algorithms. 

We  will  be  assuming  that  the  system  is  noiseless  in  that  only  transmitted  photons  are  received  and  that 
additive  noise  in  the  receiver  can  appear  only  due  to  imperfect  reception  mechanisms.  Obviously,  in  a  real 
system,  these  assumptions  would  need  to  be  examined  and  likely  relaxed. 

Radio  frequency  (RF)  signals  are  generally  modeled  with  the  viewpoint  that  they  would  have  been 
perfectly  measurable,  if  only  the  additive  noise  at  the  receiver  had  not  corrupted  the  measurement.  On  the 
other  hand,  the  quantum  nature  of  optical  signals — even  noiseless  signals — means  that  they  are,  at  their 
core,  not  perfectly  measurable.  Errors  are  made  because  of  this  quantum  measurement  inadequacy,  not 
because  of  the  corruption  of  any  additive  noise  (unless,  of  course,  noise  is  also  added  at  the  receiver,  such 
as  by  an  optical  amplifier).  The  best  receivers  are  based  on  knowledge  of  this  property. 

The  generic  communication  system  we  will  be  discussing  is  shown  in  Figure  I.  Bits  are  presented  to  the 
system  by  an  external  data  source.  We  may  decide  to  use  error-correction  coding  to  improve  power 
efficiency.  The  encoder  maps  groups  of  source  bits  into  groups  of  encoded  “channel  bits.”  These  encoded 
channel  bits  may  be  further  mapped  into  “channel  symbols,”  which  then  direct  the  operation  of  the 
modulated  optical  signal  generator.  (In  fact,  the  bits-to-symbols  mapping  may  have  been  done  at  the  input 
to  the  encoder  or  in  the  encoder  itself.) 

To  increase  transmitted  power,  there  may  be  an  optical  amplifier.  The  (lossy)  channel  is  here  assumed  to 
include  everything — transmitter  optics,  beam  spreading  losses,  channel  attenuation,  and  receiver  optics — 
between  this  amplifier  and  the  receiver,  which  then  makes  measurements  of  the  channel  symbols.  Finally, 
the  outputs  of  this  receiver  box  are  fed  to  a  de-mapper  and/or  decoder,  whose  outputs  are  estimates  of  the 
original  source  data  bits  and  which  are  sent  on  to  the  destination.  Modern  systems  may  combine  certain  of 
these  functions  in  various  ways,  but  we  will  use  this  block  diagram  with  the  understanding  that  variants 
will  have  obvious  analysis  extensions. 
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Figure  I .  The  optical  communication  system. 


The  structure  of  this  report  is  as  follows: 

Preliminaries 

•  System  Options  and  Requirements 

•  The  Classical  Receivers 

Information  Theory  of  Classical  Systems 

•  Notation  and  Basic  Error  Probability  Formulas  for  Classical  Systems 

•  Classical  Channel  Capacity 

•  Classical  Computation  Cutoff  Rate 

Information  Theory  of  Quantum  Systems 

Modulation  Formats 

•  Generalized  On-Off  Keying 

•  Orthogonal  Modulations — PPM  and  FSK 

•  Af-ary  Phase  Shift  Keying  (MPSK) 

Comparisons 

For  those  readers  interested  only  in  the  performance  results,  the  Preliminaries  sections  can  be  skipped. 
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2.  Preliminaries 

2.1  Communication  Systems  Options  and  Requirements 

Possibilities  for  signal  generation  depend  upon  the  existence  of  technologies  that  can  vary,  at  the  desired 
symbol  rate,  the  parameter  of  interest  and  then  can  send  the  constructed  signal,  at  adequately  high  power, 
to  the  receiver.  Thus,  the  selection  of  formats  will  be  directed  by  the  availability  of,  for  example,  direct 
laser  modulation  of  amplitude,  wavelength,  or  phase;  or  of  external  modulation  of  amplitude, 
polarization,  or  phase.  Further,  the  availability  of  peak-  vs  average-limited  power  generation  will  favor 
certain  choices. 

Once  a  transmission  format  has  been  selected,  there  exists  a  hierarchy  of  receiver  structures  that  can  be 
employed,  depending  on  what  is  easily  known  or  measurable  at  the  receiver.  Selection  here  will  depend 
upon  knowledge  of  signal  frequency  (or  equivalently  wavelength),  phase,  amplitude,  spatial  properties, 
waveform  shape,  timing,  and  polarization.  This  selection  is  also  driven  by  available  receiver  technologies 
with  adequate  bandwidths. 

Each  pairing  of  a  modulation  format  with  a  receiver  structure  has  its  own  set  of  analytical  tools  and 
derived  performance.  To  walk  through  this  tree  of  possibilities,  in  this  work  we  will  first  examine  the 
analysis  mathematics  applicable  to  the  various  receivers  and  then  apply  them  to  each  major  family  of 
modulations. 

Therefore,  let  us  first  survey  the  main  receiver  types. 

2.2  The  Classical  Receivers 
Noiseless  Photon  Counter 

In  this  theoretical  device,  we  make  the  classical  assumption  of  photons  being  detected  with  Poisson 
statistics  at  a  rate  proportional  to  the  instantaneous  signal  power. 

‘  =  0.1.2 . 

Here,  PPois  is  the  probability  of  counting  k  photons  in  a  measurement  period  where  N  is  the  average 
number  of  photons  (possibly  fractional)  that  could  have  been  expected  in  such  a  period.  The  perfect 
photon  counter  has  a  detection  efficiency  of  I.  (We  also  assume  in  this  work  that  there  are  no  noise 
photons  in  such  a  system.  Since  performance  in  noise  can  be  shown  to  depend  on  two  parameters — the 
signal  photons  per  symbol  and  the  noise  photons  per  symbol — and  not  just  their  ratio,  as  in  RF  systems, 
we  will  not  cover  this  case  in  this  work.  Technologies  presented  here,  however,  have  direct  applicability 
to  the  noisy  case,  with  noises  also  modeled  as  photons  counted  with  Poisson  statistics.)  Thus,  at  an 
average  arrival  rate  of  X  photons  per  second  over  a  time  of  T  seconds,  we  see  that  N  =  XT.  We  can  also 
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see  that  there  is  always  a  finite  probability  of  not  detecting  a  photon  (i.e.,  counting  0  photons)  when  a 
signal  is  actually  present,  and  this  probability  is  given  by 

P«„,(0|  N)=e-r  (2) 

(Note  that  for  a  signal  with  a  received  power  such  that  there  is,  on  average,  one  photon  to  be  detected  per 
period  of  time,  if  we  count  photons  over  many  such  periods  of  time,  we  will  count  zero  in  1/e  or  37%  of 
them.)  In  modulations  based  on  the  existence  or  non-existence  of  pulses,  then,  errors  with  such  a  receiver 
will  be  made  only  if  pulses  were  “supposed  to”  exist  but  were  not  detected.  (This  fact  is  sometimes  called 
the  quantum  limit,  but  we  will  see  that  there  are  better  uses  for  this  term.) 

True  noiseless  photon-counting  devices  with  high  detection  efficiency  and  wide  bandwidth  are  difficult  to 
build.  Technology  for  them  is  rapidly  growing,  however. 

Heterodyne 

The  architecture  for  a  (balanced)  heterodyne  front  end  is  shown  in  Figure  2. 
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Figure  2.  Heterodyne  receiver  (balanced  configuration). 


At  the  receiver,  a  strong  local  oscillator  (LO)  laser  beam  is  combined,  in  a  beam-splitter,  with  the  very 
small  incoming  optical  signal.  The  LO  is  constant  amplitude  (or  constant  wave,  CW)  and  its  wavelength 
is  selected  (or  controlled)  to  be  near  but  not  equal  to  the  wavelength  of  the  received  signal.  (Note: 
polarization  must  also  be  matched  in  classical  heterodyne  receivers,  but  there  are  variants  that  do  not 
require  polarization  [14).)  The  combined  signal  is  detected  at  (a  differenced  pair  of  balanced)  square-law 
devices  whose  electrical  output  can  be  shown  to  be  a  sum  of  the  “beat  signal”  between  the  two  signals 
plus  the  so-called  shot  noise  induced  by  the  photons  in  the  small,  received  signal.  By  the  selection  of  the 
wavelength  offset,  the  electrical  beat  signal  comes  out  at  a  desired  intermediate  frequency.  With  proper 
choice  of  this  frequency,  standard  RF  radio  receiver  techniques,  such  as  electrical  filtering,  etc.,  can  be 
employed  in  further  processing.  Thus,  this  heterodyne  front  end  can  be  followed  by  a  simple  square-law 
detector,  by  frequency-selective  elements,  or  by  coherent  receivers,  all  in  the  electrical  domain. 
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When  the  local  oscillator  power  is  very  large  compared  to  the  received  signal  power,  the  output  electrical 
signal  is  well  modeled  as  an  undistorted  but  amplified  copy  of  the  input  optical  signal,  now  embedded  in 
white  Gaussian  noise  and  with  its  carrier  frequency  now  set  at  the  intermediate  frequency.  Since  phase 


(relative  to  the  LO)  has  been  preserved,  either  coherent  or  noncoherent  techniques  can  be  used  in  further 


processing  and  in  the  detection  process. 

It  is  well  known  that  the  optimum  processing  to  perform  in  additive  white  Gaussian  noise  is  to  apply  a  so- 
called  matched  filter,  matched  to  the  signal,  with  the  output  of  the  filter  then  sampled  at  the  point  of 
largest  signal-to-noise  ratio  (SNR).  (See,  for  example,  [17].) 

This  sampled  output  has  an  SNR  equal  to  (as  described  in  RF  systems)  EJN0  (signal  symbol  energy 
divided  by  the  noise  spectral  density).  It  can  be  shown  that  for  our  heterodyne  receiver,  this  value  is  equal 
to  the  input  photons-per-symbol!  With  this  equivalence,  we  see  that  all  discussions  and  analyses  of 
classical  receivers  in  this  work  are  applicable  to  RF  systems. 

The  probability  density  function  of  the  sampled  output  signal,  properly  normalized,  at  the  output  of  a 
coherent  receiver  is  Gaussian 


fG(x-^2Ns)  , 


(3) 


where 


(4) 


and  its  cumulative  distribution  function  is 


(5) 


There  are  many  algorithms  for  efficiently  calculating  this  well-known  function. 
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A  constraint  on  heterodyne  systems  is  that  the  wide  bandwidth  electrical  lasercom  signal  must  appear  on  a 
carrier  at  an  intermediate  frequency,  which  means  that  the  detection  and  the  electrical  processing  must  be 
done  with  a  total  bandwidth  likely  1.5-3  times  as  wide  as  the  lasercom  signal  itself.  This  is  illustrated  in 
Figure  3. 

Receiver  BW 
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Figure  3.  Heterodyne  signed  structure. 


Another  constraint  is  that  the  local  oscillator  must  be  strong  enough  (a)  to  allow  the  output  beat  signal  to 
well-approximate  an  undistorted  signal  in  noise  and  (b)  to  allow  the  heterodyne  “gain"  to  be  large  enough 
that  the  Gaussian  noise  level,  set  by  the  incoming  signal,  is  larger  than  any  intrinsic  thermal  noises  in  the 
detection  and  post-processing  electronics. 

For  so-called  quadrature  modulations — meaning  that  the  signal  varies  in  the  two  dimensions,  in-phase 
and  quadrature — it  turns  out  we  will  need  to  calculate  the  two-dimensional  Gaussian  probability  of  the 
plane  missing  a  wedge-shaped  section,  as  shown  in  Figure  4. 


Figure  4.  Definition  of  Craig  fs  function. 
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Although  there  have  been  many  versions  of  formulas  defining  this  probability,  Craig’s  function  [7], 
which  can  be  expressed  as 


1  x-v 

fc  (P,V)  =  —  1  exp 

*  J 

0 

seems  to  have  the  best  properties,  including  its  ease  of  numerical  integration  and  the  fact  that  the  SNR 
value  is  presented  linearly  in  the  exponent  of  the  integrand.  Here,  ft  I 2  is  the  symbol  signal-to-noise  ratio. 

Homodyne 

The  architecture  for  a  homodyne  front  end  is  shown  in  Figure  5. 
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Figure  5.  Homodyne  receiver  (balanced  configuration). 

The  difference  between  homodyne  and  heterodyne  is  that  the  homodyne  local  oscillator  is  at  exactly  the 
same  wavelength/frequency,  polarization,  and  optical  phase  as  the  incoming  signal.  The  output  beat 
signal  is  thus  a  coherently  derived  baseband  signal  (near  DC).  There  is  no  intermediate  frequency  carrier. 
The  expanded-bandwidth  requirement  of  the  heterodyne  system  is  thereby  avoided.  Furthermore,  it  can  be 
shown  mathematically  that  the  signal-to-noise  ratio  at  the  output  of  the  ft lter/sampler  is  double  that  of 
heterodyne — that  is,  3  dB  better — due  to  the  lack  of  an  “image”  signal  at  minus  frequencies.  (See  Figure  3.) 
The  probability  density  function  of  the  output  of  a  homodyne  receiver  is  the  same  as  that  of  heterodyne,  but 
as  if  the  SNR  were  twice  as  large. 

The  drawbacks  of  homodyne  are  that  (a)  an  optical  phase-locked  loop  with  very  tight  coupling  back  to  the 
input  is  required  to  keep  the  LO  in  phase  with  the  received  signal  and  (b)  the  “baseband”  output  really 
only  carries  amplitude  and  sign  information  and  thus  supports  only  one-dimensional  signals.  Variations  of 
this  architecture  that  try  to  measure  two-dimensional,  phase-related  modulations  (such  as  QPSK — where 
a  homodyne  QPSK  receiver  splits  the  signals  with  an  optical  hybrid  [3])  turn  out  to  give  performance  that 
is  no  better  than  heterodyne.  On  the  other  hand,  such  a  homodyne  QPSK  receiver  allows  for  a  coherent 
QPSK  receiver  with  only  baseband  electronics,  a  significant  advantage  over  QPSK  heterodyne  with  its 
large  receiver  bandwidth  requirement. 
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Preamplification 


A  perfect  optical  preamplifier  is  very  like  an  RF  amplifier:  it  takes  a  very  small  (noiseless,  in  this 
discussion)  optical  signal  and  amplifies  it  with  no  distortion  up  to  a  more  usable  level,  with  the  cost  of  the 
addition  of  a  certain  amount  of  noise.  This  output  signal  level  is  defined  as  that  level  such  that  further  electro- 
optical  processing  adds  only  a  negligible  amount  more  noise.  As  in  RF  amplifiers,  noises  in  the  amplifiers  we 
will  discuss  (e.g.,  doped  fiber  amplifiers)  are  well  modeled  as  additive,  white  Gaussian  processes. 

We  will  assume  the  most  efficient  optical  amplifiers  with  3-dB  noise  figures.  This  can  be  shown  to 
correspond  to  the  fact,  as  in  the  heterodyne  receiver,  that  the  total  signal-energy-to-noise  ratio  in  a  match- 
filtered  post-amplifier  symbol  measurement  is  equivalent  to  the  average  number  of  symbol  photons  seen 
at  the  input  [22]. 

Note  that  certain  properties  of  an  optically  amplified  signal  are  the  same  as  those  of  the  electrical 
heterodyne  signal:  the  desired  signal  is  at  a  level  that  is  easier  to  detect  with  many  types  of  detectors  and 
the  otherwise  noiseless  signal  is  now  embedded  in  white  Gaussian  noise.  However,  the  heterodyne  signal 
has  been  converted  to  an  electrical  signal  on  an  IF  carrier.  The  preamplified  signal  is  still  optical  at  its 
original  wavelength.  Usually,  the  detection  process  is  noncoherent  because  if  a  coherent  method  had  been 
available,  then  a  coherent  optical  receiver  would  likely  have  been  available.  We  should  also  note  that 
optimal  matched  filtering  in  an  optically  amplified  system  must  be  in  the  optical  domain  since  only  a 
certain  fraction  of  the  noise  effects  can  be  removed  post-detection.  Usually,  it  is  difficult  to  achieve  true 
optical  matched  filtering,  although  adequate  approximations  are  available  [6]. 

We  should  make  one  important  note  here.  The  small,  noiseless,  but  inherently  unmeasurable  input  optical 
signal  has  properties  very  different  from  an  amplified  version  of  it,  which  has  had  classical  Gaussian 
noise  added  to  it.  We  will  discuss  this  further  in  Section  4. 

A  signal  that  is  preamplified,  optically  match-filtered,  envelope-detected,  and  sampled  has  the  Rician 
probability  density  function 


(7) 


where  r  is  the  envelope  of  the  signal  being  measured  (r  is  the  power)  and  where  the  envelope  of  the 
desired  signal  measurement  is 


(8) 


8 


The  randomness  exists  because  of  the  narrowband  (due  to  the  matched  filtering)  Gaussian  noise  in  the 
signal.  When  there  is  no  signal,  only  noise,  this  density  function  becomes  Rayleigh, 


f  Ray  (r)  —  f  Ric  r) 


(9) 


The  cumulative  distribution  function  of  the  Rician  density  is  called  the  Marcum  Q  function 


Q(a,b)  =  jfltk  (a,r)dr  . 


(10) 


h 


There  are  a  number  of  published  algorithms  for  efficiently  calculating  this  function  [21]. 

In  the  case  of  quadrature  modulations  with  noncoherent  detection,  we  will  need  to  calculate  the 
probability  of  measuring  phase  differences  over  wedge-shaped  areas  of  the  2-D  signal  plane.  With  many 
similarities  (by  design)  to  Craig's  function  above,  Pawula  [15]  showed  that  the  relevant  probability  can  be 
written  as 


GO 
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This  function  is  also  amenable  to  numerical  integration. 

Real  Detectors 

Options  discussed  so  far — photon  counting,  preamplified,  and  heterodyne  and  homodyne — assume  perfect, 
noiseless  detectors.  Real-world  detectors  add  noise  to  varying  degrees.  Avalanche  photodiodes  (APDs),  for 
example,  can  operate  like  photon  counters  but  with  appreciable  noise.  (Generically,  systems  that  make 
measurements  of  signal  levels,  and  possibly  add  noise,  are  known  as  Direct  Detection  systems.)  Although 
real  detector  models  should  be  used  when  analyzing  real  receivers,  their  performance  is  not  fundamental  and 
depends  on  device  design  parameters  that  trade  off  certain  inherent  noise  properties.  Therefore,  we  will 
restrict  ourselves  to  the  basic  receiver  models.  For  communications  system  variants  that  require  analyses 
with  other  noise  or  distortion  features,  it  is  left  to  the  reader  to  extrapolate  from  the  results  given  here. 

2.3  Requirements  of  the  Classical  Receivers 

We  can  see,  from  the  previous  discussions,  that  there  exists  quite  a  range  of  receiver  structures  that  can  be 
employed.  Each  receiver  type  has  certain  requirements  for  what  it  needs  to  know,  deduce,  or  control.  In 
this  section,  we  will  examine  more  closely  these  requirements. 
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The  one  universal  requirement  is  timing.  All  receivers  need  to  know,  or  deduce,  symbol  timing.  In  fact,  in 
our  noiseless  system,  a  photon  counting  receiver  can  do  without  any  other  knowledge,  although  practical 
limitations  mean  that  the  wavelength  must  be  within  the  bandwidth  of  the  detector. 

A  preamplified  receiver  requires,  in  addition,  knowledge  or  control  of,  to  some  degree,  (a)  wavelength — 
in  order  to  fall  in  the  band  of  amplification  and  the  matched  filter;  (b)  spatial  properties — because  many 
optical  amplifiers  have  single-mode  characteristics;  (c)  polarization — because  opposite-polarization 
amplifier  noise,  if  any,  needs  to  be  filtered  out  for  best  performance;  and  (d)  baseband  waveform  shape — 
because  optimal  filtering  of  the  post-amplifier  Gaussian  noise  requires  a  matched  filter  based  on  this  shape. 

Heterodyne  and  preamplified  receivers  require  knowledge  and/or  control  of  wavelength,  polarization,  and 
spatial  properties  in  order  to  generate  a  local  oscillator  that  mixes  with  the  signal.  Homodyne  receivers 
require,  in  addition,  knowledge  and/or  control  of  the  optical  phase.  As  in  the  preamplified  receiver,  the 
added  Gaussian  noise  requires  a  filter  matched  to  the  baseband  waveform,  although  in  the  heterodyne  and 
homodyne  cases,  as  mentioned  above,  the  filter  is  electrical. 

We  will  assume  throughout  our  analyses  that  fully  matched  local  oscillators  and  matched  filters  are 
available  in  our  receivers,  in  addition  to  perfect  symbol  timing. 

Because  the  homodyne  receiver  uses  the  phase  of  the  incoming  signal  and  measures  only  signal 
amplitude,  it  is  a  so-called  coherent  receiver,  to  be  analyzed  later.  Since  the  output  of  the  heterodyne 
system  is  an  electrical  IF  signal,  its  follow-on  electronics  can  employ  either  coherent  detection 
techniques — using  the  deduced  electrical  phase  of  the  IF  carrier — or  noncoherent  techniques — measuring 
only  electrical  power  using  square-law  detectors.  Once  again,  both  of  these  coherent  techniques  require 
electrical  matched  filters. 

We  must  next  decide  how  we  wish  to  make  and  process  our  measurements.  If  the  modulation  requires 
deducing  the  presence  or  absence  of  a  pulse,  there  are  several  possibilities.  In  a  noiseless,  photon¬ 
counting  system,  counting  any  photons  vs  counting  no  photons  makes  the  decision  simple.  In 
preamplified,  heterodyne,  and  homodyne  systems,  the  added  Gaussian  noise  means  that  optimum 
detection  of  a  pulse  requires  both  a  matched  filter  and  a  threshold  that  divides  outputs  into  those 
announced  as  “pulse”  and  those  announced  as  “no  pulse.”  The  optimum  threshold  must  be  based  on 
knowledge  of  the  signal  and  noise  levels,  which  both  must  be  either  measured  or  known  a  priori. 

Optimum  quantum  receivers,  to  be  discussed  at  length  later,  will  be  seen  to  require  knowledge  of  all  of 
the  parameters  discussed  so  far. 

Finally,  once  the  signal  has  been  “received,”  a  receiver  can  then  select  one  of  several  types  of 
measurements,  summarized  in  Figures  6  and  7.  If  there  is  no  decoding  to  be  done,  bits  or  symbols  must  be 
decided  upon,  in  so-called  “hard  decisions.”  If  there  is  a  decoder,  a  number  of  options  are  possible. 

For  a  decoder  to  make  maximum-likelihood  (ML)  decisions,  it  uses  all  the  information  it  gets  from  the 
receiver.  If  hard  decisions  are  all  it  receives,  then  the  decoder’s  ML  decision  is  equivalent  to  assuming  a 
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constant  average  signal-to-noise  ratio  (SNR)  through  all  the  symbols.  If  the  receiver  can  make  other 
measurements  of  the  received  waveform,  then  the  decoder’s  goal  is  to  use  this  information  in  its  ML 
decisions. 

The  optimum  receiver  calculates  the  likelihood  of  all  possible  code  words  (possibly  a  huge  number),  each 
of  which  may  last  many  bits.  The  decoding  decision  is  made  by  choosing  that  code  word  that  had  the 
greatest  (maximum)  likelihood. 

In  classical  communications  systems,  these  likelihoods  can  be  built  up  by  first  making  amplitude  and/or 
phase  measurements  independently  on  each  channel  symbol.  (These  are  sometimes  called  “soft  decisions” 
because  they  are  measurements  of  information  concerning  the  symbols,  without  actually  making  the 
“hard”  bit  decision.).  Then,  a  processor  constructs,  usually  implicitly  via  complex  mathematics,  the 
likelihood  estimates  of  all  possible  received  code  words  and  then  finally  makes  its  code  word  decoding 
decision.  Both  hard-  and  soft-decision  techniques  are  available  to  us  for  photon-counting,  heterodyne, 
homodyne,  and  preamplified  systems.  We  will  see,  however,  that  quantum  optimum  systems  have 
different  properties  in  each  of  these  steps. 


Binary 


Hard  decision 


j 


Data  bits  sent  to  destination, 
or 

Channel  bits  sent  to  decoder 


Information  sent  to  soft- 
decision  decoder 


Figure  6.  Binary  receiver  options . 


M-  ary 


Hard  decision 


j 


Soft  decision 


Measurements  of 
ail  relevant 
filter  outputs 


Symbols  sent  to  decoder, 
or 

Symbols  translated  to  bits 

•  Data  bits  sent  to  destination, 
or 

•  Channel  bits  sent  to  decoder 


Information  sent  to  soft- 
decision  decoder 


Figure  7.  M-arx  receiver  options. 


All  of  the  requirements  for  the  different  receivers  are  summarized  in  Table  I.  The  table  shows  which 
receiver  types  need  knowledge  or  control  of  which  parameters.  It  also  lists  the  subsystems  typically  used 
to  estimate  the  parameters.  The  modulations  assumed  include  OOK,  orthogonal  modulations,  and  PSK, 
all  of  which  will  be  described  in  Section  5. 
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Table  1 

Knowledge  Requirements  of  Communications  Receivers 


Feature 

Symbol  Timing 

Waveform  Shape 

Signal  Amplitude 

Receiver  Noise  Power 

Optical  Frequency/Wavelength 

Subcarrier/IF  Phase 

Optical  Phase 

Polarization  State 

Spatial  Properties 

Noiseless  Photon  Counting  Hard  Decision 

X 

Noiseless  Photon  Counting  Soft  Decision 

X 

Heterodyne  Coherent  Hard  Decision 

X 

X 

A 

A 

X 

X 

X 

X 

Heterodyne  Coherent  Soft  Decision 

X 

X 

X 

X 

X 

X 

X 

X 

Homodyne  Coherent  Hard  Decision 

X 

X 

A 

A 

X 

X 

X 

X 

Homodyne  Coherent  Soft  Decision 

X 

X 

X 

X 

X 

X 

X 

X 

Heterodyne/Pre-Amp  Non-Coherent  Hard  Decision 

X 

X 

A 

A 

X 

X 

X 

Heterodyne/Pre-Amp  Non-Coherent  Soft  Decision 

Lx 

X 

X 

X 

X 

X 

X 

(Noiseless)  Quantum  Hard  Decision 

X 

X 

X 

X 

X 

X 

X 

(Noiseless)  Quantum  Max  Likelihood 

X 

X 

X 

X 

X 

X 

X 

Means  to  obtain  knowledge 

(A  =  OOK  only) 

Time  tracking  loop 

a  priori  knowledge 

Power  monitor 

Out-of-band  monitor;  Demodulator/decoder  outputs 

Freq-tracking  loop;  a  priori  via  calibration 

IF  phase  lock  loop 

Optical  phase  lock  loop 

Polarization  tracker/controller 

WF  sensor;  camera 

Now  that  we  have  finished  this  whirlwind  tour  of  receivers  and  have  introduced  the  mathematical 
functions  needed  for  their  analyses,  we  will  examine  the  metrics  of  communication  system  performance. 
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3.  Information  Theory  of  Classical  Systems 

3.1  Notation  and  B  asic  Error  Probability  Formulas  for  Classical  Systems 

We  will  need  to  distinguish  between  peak  and  average  powers  as  well  as  energy  per  pulse,  bit,  and 
symbol,  so  we  propose  the  following  notation  throughout:  Whether  a  channel  symbol  corresponds  to  a 
single  bit  or  multiple  bits,  we  will  call  it  “channel  symbol.”  Let  NP  be  the  average  number  of  photons  in  a 
received  pulse,  Ns  the  average  number  in  a  channel  symbol,  and  Nh  the  average  number  in  the  time 
corresponding  to  a  source  bit.  Let  TP  be  the  time  duration  of  a  pulse,  7$  the  time  duration  of  a  channel 
symbol,  and  Th  the  time  duration  of  a  source  bit.  We  thus  have  that  NP/TP  is  the  (pulse-averaged)  peak 
received  power  and  N$/Ts  is  the  average  received  power.  Relationships  between  the  various  parameters 
will  be  explained  for  each  format. 

For  soft-decision  systems,  we  use  the  notation  PqxO’l*)  to  represent  the  probability  density  function  of 
receiving  the  (sampled)  random  quantity  Y  =  y  given  that  the  Ath  symbol  at  the  input  X  was  sent.  (We  can 
also  see  that  an  even  more  general  case  allows  nonquantized  inputs,  X).  We  may  drop  the  subscripts  X  and 
Y  when  the  variables  are  obvious.  The  value  y  may  be  the  single  voltage  sample  at  the  output  of  a 
matched  filter,  or  it  may  correspond  to  a  vector  of  measurements  of  M  sampled  filter  outputs.  We  number 
these  symbols  from  0  to  M  -  I .  For  hard  decision  systems,  we  use  the  notation  Py\x(i\k)  to  represent  the 
probability  that  symbol  i  is  announced  by  the  receiver  when  symbol  k  is  sent.  For  M- ary  and  coded 
systems,  we  use  the  notation  p,  to  mean  the  input  probability  (or  frequency,  or  prior  probability)  that 
symbol  i  is  sent. 

The  average  probability  of  a  hard  decision  M-ary  symbol  error  is  given  by  the  formula 

m- i  m- i  (12) 

ps  =  Z  2>*ph*oi*) 

7=0  k=0.* j 

=  ]~Y^Pjpy\xU\j), 


which,  for  binary  systems,  is  equivalent  to 

Ps=PlP(0\\)+PoP(\\0)  ■  (13) 

Note  that,  if  p()  =  P\,  then  a  symbol  error  is  equivalent  to  a  bit  error,  Ph ,  where  we  let  Ps  and  P h  be, 
respectively,  the  symbol  and  bit-error  probabilities. 


15 


Remember  from  basic  probability  theory  that 


M 


(14) 


^  F*y\x  0  I  * 


and 


M 


(I?) 


( y I ^ )Pk  ~  Py (j ) * 


*=1 


where  we  let  P>(/)  define  the  unconditional  probability  density  function  of  the  output,  Y  =  j. 

In  Section  5,  we  will  investigate  two  basic  classes  of  modulations:  binary,  and  nonbinary.  The  nonbinary 
modulations  we  will  investigate  will  all  have  a  certain  symmetry  in  that  each  of  the  M  signals  is  related  to 
the  set  of  remaining  signals  in  the  same  way.  That  is. 


P(j\k)=  P(j+\\k  +  \) 
=  P{j-k  |0)  , 


(16) 


where  the  sum  and  difference  are  taken  modulo  M.  Since  no  signal's  relative  properties  are  different  from 
any  other’s,  we  will  be  assuming  that  the  prior  (i.e.,  the  input)  probabilities  are  equivalent,  i.e.,  pk  =  I/M 
for  all  k  in  the  symmetric  M-ary  system.  In  this  case.  Equation  (12)  would  reduce  to 


(17) 


Ps  =  1-/>(0|0). 


Furthermore,  Equation  (15)  would  reduce  to 


(18) 


which  means  that,  in  such  symmetric  systems,  the  outputs  are  also  equiprobable. 

One  last  topic  regarding  M-ary  symbols  is  the  mapping  of  bits  into  symbols  and  the  symmetrical  problem 
of  mapping  received  symbol  decisions  into  source  bits.  It  is  of  interest  to  know  the  error  probability  of  the 
final  bits  when  it  was  the  symbol  error  probabilities  that  are  usually  calculated.  Of  course,  every  system’s 
unique  mapping  needs  to  be  examined  on  its  own  merits.  However,  in  many  mappings — in  particular  so- 
called  Gray-code  mapping,  where  there  is  exactly  one  bit  difference  between  subsequent  symbols — it  can 
be  shown  that  the  probability  of  a  bit  error  is 


(19) 
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3.2  Classical  Channel  Capacity 


Definition 

Channel  capacity  is  the  theoretical  maximum  rate  at  which  source  data  can  be  transmitted  error-free  over 
a  noisy  or  quantum  channel  to  the  destination.  It  is  the  modern  metric  against  which  all  communications 
system  performance  should  he  compared.  Capacity,  C,  is  a  function  of  the  average  number  of  photons 
per  channel  symbol  that  is  available  at  the  destination.  It  is  usually  calculated  in  (source)  bits  per  channel 
symbol  time,  a  rate,  but  can  also  be  specified  in  source  bits  per  second  or  some  related,  correctly 
normalized  quantities. 

It  is  known  that  certain  long,  complex  coding  schemes  exist  that  can  allow  a  communication  system  to 
approach  this  performance  limit.  A  code  at  rate  R  <  I  means  that  the  code  transmits  I  channel  bit  for 
every  R  source  bits.  If  we  employ  a  capacity-achieving  code,  i.e.,  a  highly  efficient  code  at  code  rate,  C, 
then,  by  definition,  it  will  be  able  to  achieve  error-free  transmission  through  the  channel  at  the  theoretical 
minimum  power  per  original  source  bit. 


The  error  probability  of  a  coded  system  is  lower-bounded  by  the  following  formula  [17], 


Ph>H 


1- 


CjRNj, 

R 


(20) 


where  the  binary  entropy  function  is  defined  as, 

H(p)  =  ~[p\og(p)  +  (]  ~/>)log(  |  -  p)]. 


(21) 


We  will  use  only  base-2  logarithms  throughout.  This  entropy  function  is  plotted  in  Chart  I 

In  Equation  (20),  C(NS)  is  the  capacity  in  source  bits  per  channel  symbol  as  a  function  of  photons  per 
symbol,  Ns,  and  R  is  the  code  rate  employed.  (Thus,  R  /V,,  =  Ns  ).  Since  H{p)  is  positive,  we  see  that  the 
error  rate  cannot  go  to  zero  if  R  >  C. 

Capacity  Formulas — Hard  Decision 

In  this  section,  we  will  give  the  basic  formulas  for  classical  channel  capacity.  We  will  present  the 
formulas  for  both  hard-  and  soft-decision  systems.  We  will  also  provide  some  special  cases,  which  we 
will  refer  to  later. 


Classical  channel  capacity  for  hard-decision  systems  can  be  calculated  as  [I  I  ] 


Cc 


max 

{/>*} 


Y^P'Py\xU\i)  log 

•j 


Pv\x  U  I ') 

Py(j) 


(22) 
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(23) 


py(j)  =  'ZdPtPv\x(j\k), 
k 


where  Py(j)  is  the  average  probability  that  output  j  is  produced. 

(We  are  using  the  term  “classical”  here  to  distinguish  it  from  the  optical  quantum  systems,  to  be  described 
in  the  next  section).  This  value  has  units  of  source  bits  per  channel  symbol.  The  formula  is  shown  as  a 
maximum  over  all  choices  of  prior  probabilities  since  this  is  usually  a  free  parameter  for  the  system 
designer. 

The  formula  can  also  be  written  in  a  way  that  shows  the  symmetry  between  the  input  and  the  output. 
(Capacity  is  based  on  the  average  mutual  information  between  the  inputs  and  the  outputs.)  The  following 
equation  is  equivalent  to  Equation  (22): 


max  log  xr  J 

Pk  i.j  W^(i) 

(24) 

f*XY  y)  =  F*Y\X  0  1 

(25) 

PX  (0  =  Pi 

(26) 

For  a  binary  input/output  channel,  we  can  rearrange  Equation  (22)  to  be 

Cr  =  max{tf  />,/>(!  |l)  +  A>0P(l|0)]-ft,//  P(l|0 )\-pxH  mil)!}. 


(27) 


In  the  case  of  symmetric  M-ary  modulations,  the  formula  for  channel  capacity  reduces  to 

Cc  =  log(M)  +  X/5(;'|0)log/>(y|0), 

j 


(28) 


which,  when  describing  the  binary  symmetric  channel  (BSC;  M  -  2,  p\  =  p0  =  1/2)  further  reduces  to 

Cc=l-//[/>(l|0)].  (29) 

Channel  capacity  for  the  BSC  as  a  function  of  symbol  error  probability,  P  (0|  I),  is  plotted  in  Chart  2. 
This  curve  shows  us  the  symbol  error  probability  at  which  a  channel  capacity-achieving  code  will 
“operate.”  Notice  that  the  lowest  rate  capacity-achieving  codes  can  operate  with  quite  high  crossover 
probabilities. 
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Capacity  Formulas — Soft  Decision 

For  a  soft-decision  system,  classical  channel  capacity  is  calculated  as  (with  the  most  general,  analog 
transmitted  inputs,  X) 


f  Py\ a'(v|*) 

Cc=  sup  pxU)Pm(y  I*)  log - 7 —dxdy 

(P.vU> I  J  Py(v) 

Py  (>’)  =  \Pr\x(y\x)Px(x)dx 

or,  for  a  more  typical  M- ary  input  alphabet, 

f  /  i  •  \  i  Pv\  x  i y  I '  ^  i 

Cc  =  max  2_  I  P,Py\  x  ( >’  0  log - — -  ay 

(ft)  ;  J  pY(y) 

pr(y)  =  Y,PiPY\x(y\iy 


(30) 


(31) 


(32) 

(33) 


Perhaps  a  more  useful  way  to  write  Equation  (32)  is 


c,.  -  max{-2;Jwm(.vl01ogZP, 


From  this,  the  special  case  of  binary  systems  reduces  to 


Po  +  P\ 


Cc  =  -Po  Jp(.y|0)log 

~P\  jp(.v|l)log  Pi  +  Po 


piy  1 1) 

p  (y  1 0) 

P  ( v  1 0) 

Piy  ID 


dy 


dv 


(34) 


(35) 


or,  equivalently, 

Cc  =  p,  Jp(v|  l)logp(y|  1  )dy  +  p0  Jp(.v|0)logp(v|0)</v  -  Jp(.v)  log  piy)  dy , 


(36) 


where 

piy)  =  P\Piy  ID +PoPiy  1 0). 


(37) 
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The  soft-decision  capacity  using  A/-ary  symmetric  signals  reduces  to 


cc  = 


log (M)-  |p(  y\  0)log 


1  + 


y1'  /?(  yl  A' ) 
hp(y  |0) 


(38) 


We  should  note  that  the  measurement,  y,  here  may  actually  correspond  to  a  vector  of  measurements  (such 
as  the  outputs  of  several  different  filters),  and  thus  the  integration  over  y  may  be  a  multiple  integration. 

Except  in  very  special  cases,  we  will  find  that  soft-decision  channel  capacity  needs  to  be  calculated  via 
numerical  integration  because  of  the  log  inside  the  integral.  For  the  even  more  complex  multiple-integral 
version,  a  technique  often  used  is  to  notice  that  the  integral  is  a  kind  of  expected  value.  Thus,  a  Monte 
Carlo  simulation  using  computer-generated  random  variables  can  be  used. 

It  can  be  shown  (and  was  shown  by  Shannon  in  his  seminal  work)  that  starting  with  Equation  (30)  and  in 
an  otherwise  unconstrained  additive  Gaussian  noise  channel,  the  channel  capacity  is  [18,  19] 


Cc  -  log(  I  +  NS  ). 


(39) 


(For  the  receivers  we  will  examine,  this  formula  will  be  relevant  to  the  coherent  heterodyne  architecture.) 

As  discussed  in  an  earlier  section,  the  noise  in  a  homodyne  receiver  is  also  Gaussian  but  at  half  the  level. 
Since  homodyne  systems  do  not  allow  for  quadrature  information,  however,  they  are  only  half  as 
bandwidth  efficient  as  heterodyne. 

The  capacity  curves  for  heterodyne  and  homodyne  unconstrained  signals  are  plotted  in  Chart  3. 
(Sometimes  the  heterodyne  equation  is  shown  with  other  factors  of  2,  when  the  normalization  is  with 
respect  to  dimensions.  Equation  (39)  is  correct  for  source  bits  per  channel  usage.)  This  ultimate  capacity 
can  be  achieved  using  Gaussian-distributed,  i.e.,  not  quantized  to  M  levels,  source  signal  values. 

Measure  of  Efficiency 

Although  the  basic  value  of  channel  capacity — source  bits  per  channel  symbol  at  a  given  symbol  photon 
count — is  highly  useful,  we  sometimes  want  to  know  the  cost  of  achieving  a  certain  data  rate  as  opposed 
to  knowing  the  data  rate  achievable  at  a  certain  cost  (signal  power).  A  way  to  do  this  is  to  calculate  the 
photons  required  per  source  bit  as  a  function  of  some  speed  or  redundancy  property  of  the  code.  We  will 
present  this  efficiency  at  channel  capacity  vs  bandwidth  expansion,  where  bandwidth  expansion  is  a 
combination  of  the  expansion  due  to  the  code  parity  check  overhead  plus  the  expansion  due  to 
modulation. 

Since  Ns  =  average  photons  per  channel  symbol,  and  Cc  =  source  bits  per  channel  symbol,  we  see  that 
Ns  /  Cc  =  average  photons  per  source  bit.  That  is,  Ns/Cc  is  the  number  of  photons  required  at  the  receiver 
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in  order  for  a  system  to  reliably  deliver  source  bits  using  a  channel-capacity-achieving  code  of  rate  Cc. 
The  bandwidth  expansion  due  to  coding  is  just  \/Cc,  and  the  bandwidth  expansion  from  modulation  must 
be  calculated  on  a  case-by-case  basis.  As  an  example  of  efficiency  at  channel  capacity,  see  Chart  4,  in 
which  we  have  plotted  the  efficiency  of  the  Gaussian  channel  capacity  for  both  heterodyne  and  homodyne 
systems.  Notice  that  heterodyne  efficiency  bottoms  out  at  In  2,  or -1.59  dB  (with  homodyne  3  dB  lower). 
Thus,  with  a  (high  bandwidth,  low  code  rate)  capacity-achieving  code,  we  should  be  able  to  transmit 
source  bits  at  rate  /?,,  by  delivering  to  the  heterodyne  receiver  Rh  In  2  photons  per  second.  We  can  also  see 
that  heterodyne  becomes  more  efficient  than  homodyne  for  normalized  bandwidths  less  than  about  0.63 1 . 

3.3  Classical  Computation  Cutoff  Rate 
Definition 

For  many  years,  the  computation  cutoff  rate  (sometimes  called  computational  cutoff  rate  or  even  just 
cutoff  rate)  was  considered  the  maximum  rate  achievable  by  buildable  codes.  With  the  advent  of  turbo- 
codes  and  their  progeny,  classical  channel  capacity  (which  is  greater  than  cutoff  rate)  is  now  known  to  be 
essentially  achievable,  although  the  cutoff  rate,  Rc  (sometimes  called  R0)  is  still  a  good  upper  bound  for 
many  classical  coding  schemes.  It  is  similarly  parameterized  as  source  bits  per  channel  symbol. 

Cutoff  Rate  Formulas — Hard  Decision 

The  classical  computation  cutoff  rate  for  hard-decision  systems  is  defined  by  the  formula  [II] 


2  ' 

Rc  =  —  log" 

min  V 

{pk  \  “ 

^  Pk  0  1 

k 

> 

(40) 


For  a  general  binary  channel,  this  formula  is 


l-2PiPo 


I  -  I  0)(1  -  P(0 1  I)) 
-VP(0|1)(1-P(1|0)) 


(41) 


For  our  symmetric  M- ary  channel,  we  find 


Rc  =  log(A/)  —  21og 


XVmjo) 


which,  for  the  BSC,  further  reduces  to 


Rc  =  1-log 


1+2A//,(1|0)P(0|0) 


The  cutoff  rate  for  the  BSC  is  shown  in  Chart  5. 


(42) 


(43) 
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Cutoff  Rate  Formulas — Soft  Decision 


For  a  soft-decision  decoder,  the  classical  computation  cutoff  rate  is 


Rc  =  -  log <  inf  |  [  jpx  (x)pY\x  ( .V I x  )u2dx 


cly 


for  the  most  general  analog  input  case.  For  the  more  typical  discrete,  M-ary  input  case,  it  is 


Rc  =  -  log 


min  f 

,1  J 


lPk 


<j\ 


(44) 


(45) 


For  a  binary  input  system,  this  reduces  to 

/?c=-l°g{  ]~2PoP\  +2PoP\ \yfpi y|0)p(v|  l)</y}  . 


(46) 


As  in  the  capacity  calculation,  we  should  note  that  the  measurement,  y,  may  be  a  vector  of  measurements 
and  thus  the  integration  over  y  may  actually  be  a  multiple  integration. 

We  can  see  that  the  soft-decision  cutoff  rate  is  made  up  of  the  log  of  sums  of  integrals,  each  integral  being 
the  geometric  mean  of  two  density  functions.  Such  an  integral  is  called  a  Bhattacharya  coefficient.  Let  us 
define  this  integral,  with  a  nonstandard  but  easily  understood  notation,  as 

®(f(y)fg(y))=  \>lf(y)g(y)dy  ■  (47) 
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In  the  M- ary  symmetric  case,  the  cutoff  rate  can  be  reduced  to 


Rc  =  log (M)-  log 


E«[p(y|0),  p(y\k)] 


(48) 


For  the  BSC,  the  soft-decision  cutoff  rate  further  reduces  to 

Rc  =  1 — log  [l+«[p(  v|0),/?( y|  1)]J . 


(49) 


It  can  be  shown,  starting  with  Equation  (44),  that  in  an  otherwise  unconstrained  additive  Gaussian  noise 
channel,  the  cutoff  rate  for  our  heterodyne  system  is  [20] 


R 


0 


I 

ln(2) 


+  log. 


2 


(50) 


—  R0{NS ) . 

As  in  our  capacity  calculations,  the  version  for  homodyne  is  3  dB  better,  but  half  as  bandwidth  efficient. 
These  formulas  are  plotted  in  Chart  6.  The  efficiencies  at  these  Gaussian  channel  cutoff  rates  are  plotted 
in  Chart  7.  We  see  that  they  bottom  out  at  2  In  2,  or  1.42  dB,  and  In  2  or  -1.59  dB,  exactly  3  dB  higher 
than  the  efficiencies  at  Gaussian  channel  capacity  for  heterodyne  and  homodyne. 

Calculation  of  these  cutoff  formulas  is  much  simpler  than  those  of  the  soft-decision  channel  capacity. 
Very  often,  Bhattacharya  coefficients  even  have  a  closed  form. 
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4.  Information  Theory  of  Quantum  Systems 

It  has  been  known  since  the  1960s  that  the  optimum  approach  to  optical  signal  detection  is  based  on 
quantum  state  descriptions  and  projection  operators.  As  we  are  interested  only  in  optical  systems  in  which 
there  is  a  large  loss  between  the  transmitter  and  receiver,  we  will  be  considering  only  so-called  pure, 
coherent  states.  In  these  cases,  solutions  to  the  mathematical  optimizations  have  been  found  for  most 
relevant  formats. 

Interpretation  of  the  mathematically  optimum  receivers  as  implementahle  physical  systems  is  still 
largely  unknown  and  the  topic  of  much  research . 

The  following  section  is  an  extremely  cursory  introduction  to  the  relevant  theory  and  mathematics.  (See 
[12]  or[  1 3]  for  more.) 

Definitions 

Quantum  states  are  described  using  infinite-dimensional  vectors  |  y/}  in  a  Hilbert  space  over  the  field  of 
complex  numbers,  with  Hermitian  conjugate  vectors  {if/ 1,  inner  products  (lf/\lfj  ,  outer  products 
and  matrix  trace  Tr(.)  defined  in  the  obvious  ways.  It  is  the  outer  product  that  describes  a  so-called  pure 
state. 

A  coherent  state  has  the  further  description 

00  r/"  (51) 

|a)  =exp(-|«l72 )  X  77777I") 

n=()  \n-) 

for  some  complex  number,  a,  where  |/?)are  orthonormal,  so-called  number  eigenstates.  (It  is 
straightforward  to  deduce  the  Poisson  counting  statistics  from  this  model.)  Thus,  the  inner  product 
between  two  coherent  states  can  be  calculated  to  be 

|(«|/5}|!=exp(-|a-/?|2)  (52) 

and  so 

(a\a)=\.  (53) 

Notice  that,  by  Equation  (52),  there  cannot  be  truly  orthogonal  coherent  states. 

As  in  the  classical  description,  transmitters  can  select  a  coherent  state  signal  to  transmit  with  prior 
probability  p,. 
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Measurements  at  the  receiver  are  described  by  so-called  probability  operator  measures,  which  are  sets  of 
nonnegative-definite  Hermitian  operators,  summing  to  the  identity.  For  the  optical  systems  we  are 
interested  in,  the  signal  states  are  coherent,  pure,  and  linearly  independent  [12].  In  this  case,  the  operators 
can  be  shown  to  take  the  form  of  outer  products,  called  projection  operators,  of  orthonormal  vectors, 

|  !jk } ,  which  sum  to  the  identity 


(54) 


It  is  the  task  of  the  receiver  designer  to  find  the  set  of  \r/kj  to  optimize  the  relevant  metric. 
Unfortunately,  measurements  described  by  operators  do  not  (yet)  have  simple  physical  interpretations. 
Thus,  most  of  the  results  in  this  quantum  section  are  of  theoretical  interest  only  although  the  excellent 
performance  they  predict  will  provide  a  great  incentive  to  find  such  physical  interpretations. 

As  discussed  in  Section  1.0,  noiseless  quantum  systems  incur  errors  not  because  of  intrinsic  added  noise 
but  because  the  signals  themselves  are  not  completely  measurable.  Furthermore,  a  particular  physical 
measurement  on  a  quantum  signal  means  that  the  signal  is  then  no  longer  available  for  making  further 
arbitrary  measurements.  This  is  the  major  difference  between  quantum  systems  and  RF  systems.  In  an  RF 
system,  by  including  a  low-noise  front-end  amplifier,  the  RF  signal  plus  noise  can  be  made  large  enough 
so  that  it  can  be  split  into  many  replicas,  each  of  which  can  be  measured  by  a  different  matched  filter. 

We  could  try  the  same  approach  with  optical  signals,  i.e.,  put  the  optical  signal  through  an  optical 
amplifier  so  that  we  could  make  multiple  measurements  at  its  output.  Unfortunately,  as  we  mentioned 
earlier,  this  signal  in  Gaussian  noise  would  have  substantially  different  properties  from  the  original 
noiseless  signal.  Optimum  processing  requires  a  very  different  approach  for  multiple  measurements  that 
must  be  based  on  operator  theory. 

We  can  see  that  the  fact  that  the  operators  sum  to  the  identity  specifies  how  any  set  of  multiple 
measurements  must  be  related.  Changing  one  operator  element  (“matched”  to  a  zero  signal,  say)  means 
that  the  other  elements  (“matched”  to  the  other  signals)  must  also  change.  Remember  that  in  an  RF 
system,  each  of  the  post-preamplifier  filters  may  be  selected  independently. 

Quantum  Optimum — Hard  Decision 

Hard  symbol  decisions  can  be  made  with  these  quantum  operators.  We  choose  them  to  minimize  the 
probability  of  symbol  error,  defined  as  in  Equation  (12). 

Although  the  quantum  mechanical  mathematics  is  infinite -dimensional,  it  can  be  shown  that  the  coherent 
problem  can  be  reduced  to  one  in  many  fewer  dimensions.  We  only  need  to  consider  the  variables 

xjk  =  (J1j\¥k)  <55> 


for  the  pure  state. 
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Pt=WkXvA  k  =  0,\,2,....,M  - 1 


(56) 


with  prior  probabilities  pk ,  and  the  operator 


nj=\tjj)(nJ\  j= 0,1,2, . l- i 


(57) 


with  the  property 


L-l 

ln,  =  .. 


L-l 


(58) 


7=0 


The  basic  property  of  projection  operators  is  that 


P(j\k)=\x 


(59) 


1(^7; Vl)  I2 


That  is,  the  receiver  announces  signal  j  when  k  is  true  with  probability  P(j\k)  defined  here. 

Notice  that  the  dimensionality  of  the  operators  does  not  need  to  be  the  same  as  the  dimensionality  of  the 
signals.  If  L  =  M,  then  it  is  a  true  hard  decision,  with  announcements  of  the  signal  decision.  However,  we 
can  choose  L>  M  under  some  soft-decision-like  metric,  if  we  like.  (Analogously,  in  classical  systems,  we 
can  perform  hard  decisions,  or  describe  the  outputs,  y,  with  2,  3,  or  more  bits.  The  more  bits  we  choose, 
the  closer  we  get  to  soft  decisions.)  For  our  quantum  hard-decision  calculations,  we  will  use  only  L  =  M. 
We  will  discuss  the  L>  M  case  in  the  next  section. 

It  has  been  shown  (See  1 12]  for  a  summary.)  that  to  minimize  the  probability  of  symbol  error,  the  solution 
must  satisfy  the  following  set  of  nonlinear  equations: 


(60) 


PjXkjXjj  =Pkxkkxjk  , 


where  it  is  simple  also  to  find  that 


(61) 


Even  these  equations  are  difficult  to  solve  in  general  in  closed  form  although  solutions  have  been  found 
for  our  binary  and  M-ary  symmetric  cases. 
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These  hard  symbol  decisions  can  be  used  as  is,  or  can  be  fed  to  a  classical  decoder  with  performance 
described  using  the  formulas  of  classical  channel  capacity  Equation  (22)  or  classical  computation  cutoff 
rate  Equation  (40)  with  P(/'/A)  from  Equation  (59).  This  approach  is  sometimes  called  “semiclassical” 
because  the  hard  symbol  decisions  are  made  using  optimum  quantum  techniques,  but  the  decoding  is 
done  classically. 

Because  of  the  property  mentioned  at  the  end  of  the  previous  section,  the  operators  that  optimize  the 
probability  of  error  are  not  necessarily  the  same  as  those  that  optimize  semiclassical  channel  capacity  or 
semiclassical  computation  cutoff  rate.  In  other  words,  the  optimum  hard  symbol  decisions  derived  above 
for  minimization  of  the  symbol  error  probability  may  not  be  the  same  decisions  that  optimize  capacity  or 
cutoff. 

It  turns  out  that,  for  binary  (and  Af-ary  symmetric)  systems,  it  is  known  that  the  two  (or  M)  measurements 
that  make  up  the  /Voptimizing  solutions  do,  in  fact,  satisfy  necessary  conditions  for  semiclassical  capacity 
optimization.  (As  of  this  work’s  writing,  sufficient  conditions  are  not  known.)  Cutoff  optimization  has 
similar  properties.  It  has  been  shown,  however,  that  the  optimum  capacity-achieving  solution  has  d  <L<d2 , 
where  d  is  the  dimension  of  the  Hilbert  space  in  which  the  signals  can  be  defined  [8]. 

Quantum  Capacity — Maximum  Likelihood  (Soft  Decisions) 

As  discussed  in  the  previous  section,  we  can  choose  L,  the  number  of  operators  in  Equation  (57)  to  be 
greater  than  Af,  the  number  of  symbols.  This  is  analogous  to  soft  decisions  in  RF  systems,  where  each  of 
M  measurements  might  have  2  or  more  bits  of  resolution  in  them,  resulting  in  more  than  A/  possible 
outputs.  All  these  bits  of  descriptions  can  then  be  fed  to  a  soft-decision  decoder.  We  will  see  that  some 
A/-ary  quantum  receivers  do,  indeed,  perform  better  in  coded  systems  when  L>  M. 

In  the  correct  analogy  with  the  classical  maximum-likelihood  decoder,  though,  the  fully  optimum 
quantum  decoder  tries  to  make  ML  decisions  between  all  the  possible  received  codewords.  As  an 
example,  for  an  (/7,  k)  block  code  with  2k  codewords,  each  constructed  of  n  >  k  binary  channel  symbols, 
the  quantum  decoder  views  the  problem  as  an  M  =  2*-ary  problem.  As  mentioned  earlier,  in  classical 
systems,  this  is  done  implicitly  by  measuring  each  channel  symbol  and  then  computing  the  2k  ML 
estimates  after  the  fact  in  a  processor.  In  quantum  systems,  however,  this  estimate  cannot  be  built  up  by 
measuring  each  sequential  symbol  independently,  but  must  be  done  by  making  measurements  on  the 
complete  codeword  as  if  it  were  a  much  more  complex  quantum  state.  This  is  sometimes  known  as  an 
‘'entangled”  quantum  receiver. 

It  will  be  seen  that  such  an  optimum  quantum  receiver  has  a  strictly  better  performance  than  any  other 
receiver,  including  the  optimized  semiclassical  decoding  receiver.  Its  performance  is  described  thus: 
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Calculate  the  following  description  (density  operator)  of  the  quantum  information  source 

M 

Pp  =  I>*  WkYWk  I  • 

*=i 


(62) 


where  we  use  the  subscript  p  to  remind  us  that  it  is  a  function  of  the  prior  probabilities,  which  can  be 
selected  to  optimize  performance.  The  quantum  channel  capacity  can  be  shown  to  take  the  form  [13] 


Cn  =  max  —Tr 

y  (ft) 


Pp  log Pp 


(63) 


(which  is  known  as  the  von  Neumann  entropy).  It  is  usually  calculated  by  first  finding  the  eigenvalues,  A„ 
of  p  .  By  Equation  (62)  we  see  that  these  are  nonnegative  and  sum  to  I,  and  are  thus  a  lot  like 
probabilities.  Using  these  eigenvalues,  we  can  equivalently  calculate 

u  (64 

ce  =  -£a4  iog/t*  , 

k= 1 


where  we  have  suppressed  the  maximization  over  the  priors. 

For  Af-ary  modulations  with  the  symmetry  property,  we  will  find  it  useful  to  use  the  related  quantity 

vk  =  MAk ,  (65) 


which  leads  to 


M 


CQ  =logM-  —  logv 


(66) 


It  can  also  be  shown  that,  for  this  M- ary  symmetric  case,  the  symbol  error  probability  described  by 
Equation  (60)  can  be  written  explicitly  as 


Ps  = 


(67) 


Similar  to  the  classical  case  of  unconstrained  signals  in  Gaussian  noise,  there  is  a  so-called  “Gaussian" 
quantum  capacity  for  unconstrained  coherent  signals,  given  by  [13]: 

CQ  =  (\  +  Ns)\og(\  +  Ns)-Ns\ogNs  -  (68) 
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This  formula  is  plotted  in  Chart  8.  It  is  the  ultimate  performance  of  all  the  methods  we  will  have 
examined. 

We  have  also  plotted  the  efficiency  at  Gaussian  quantum  capacity  in  Chart  9.  The  most  interesting  feature 
of  this  curve  is  that  it  does  not  bottom  out  as  does  efficiency  at  Gaussian  classical  capacity.  Instead,  it 
continues  to  improve  as  more  and  more  bandwidth  is  used — that  is,  as  the  error-correction  coding  rate 
gets  smaller  and  smaller.  Theoretically,  the  channel  capacity  can  be  shown  to  become  infinite ,  i.e.,  the 
efficiency  approaches  an  infinite  number  of  bits  per  photon,  although,  of  course,  real-world  (and  perhaps 
even  model  validity)  constraints  prevent  this.  Practical  limits  on  efficiency  seem  to  be  about  10,  or 
perhaps  20,  bits  per  photon. 

The  fact  that  channel  capacity  (in  source  bits  per  channel  symbol)  increases  as  the  receiver  entangling 
increases  from  1  to  2  to  n  consecutive  channel  symbols  is  sometimes  called  superadditivity. 

Quantum  Cutoff  Rate  -  Maximum  Likelihood  (Soft  Decisions) 

The  quantum  computation  cutoff  rate  is  described  similarly.  Using  the  same  source  description  as  in 
Equation  (62),  it  can  be  shown  that  [1,  13] 


Rq  S  -  max  logC Trp  ] 

(/>J 

M 

=  -  max  log  2 

(/’<)  ft 


For  symmetric  modulations,  we  have 


rq  =  log  M  ~  log  ~~ 


(69) 


(70) 


Similar  to  the  classical  case  of  unconstrained  signals  in  Gaussian  noise,  there  is  a  “Gaussian”  quantum 
cutoff  rate,  given  by  [  13]: 


R„  = 


ln(2) 


1  +  2 Ns  ~-\[\  +  4NS2 


+  log. 


i(l+Vl  +  4N,! 


(71) 


=  R0(4NS )  . 


(72) 


[Compare  Equation  (50).] 
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The  Gaussian  quantum  cutoff  rate  is  plotted  in  Chart  10,  and  the  efficiency  at  Gaussian  quantum  cutoff 
rate  is  plotted  in  Chart  1 1 .  We  see  that  this  efficiency  is  exactly  6  dB  lower  than  that  of  Gaussian 
(heterodyne)  cutoff,  and  thus  bottoms  out  at  -4.6  dB.  Thus,  there  is  a  big  incentive  (more  than  5  dB)  to 
find  coding  schemes  that  achieve  quantum  capacity  and  not  just  quantum  cutoff  rate. 

The  Square-Root  Measurement 

Although  Equations  (60)  and  (61)  are  difficult  to  solve  in  general,  a  set  of  solutions  has  been  proposed 
that  satisfies  these  equations  in  many  cases.  Let  us  define  measurements  as 


where  the  matrix  inverse  may  need  to  be  a  generalized  inverse.  This  has  come  to  be  known  as  the  square- 
root  measurement  and  has  been  investigated  in  several  works  as  good  approximations  to  optimum 
solutions.  In  fact,  it  was  recently  shown  [10]  to  be  symbol-error-rate  optimum  for  a  wide  class  of 
modulations,  including  all  binary  modulations,  as  well  as  equal-priors  orthogonal  and  MPSK. 

Conclusions 

We  end  this  section  by  reminding  the  reader  once  again  that  physical  implementations  of  quantum 
optimum  receivers,  demodulators,  and  decoders  are  still  unknown,  as  are  codes  that  might  use  them. 
These  are  all  topics  of  much  research. 

(It  would  be  fair  to  note  that  some  years  ago  Sam  Dolinar  at  MIT  invented  a  nearly  implementable 
receiver  that  achieves  the  optimum  quantum  bit-error  probability  for  any  binary  modulation.  (See  [12]  for 
a  simplified  description.)  This  receiver  has  certain  similarities  to  homodyne  systems  in  that  it  requires  a 
phase-coherent  LO.  Unfortunately,  neither  Dolinar  nor  subsequent  researchers  have  successfully 
generalized  this  method  in  an  implementable  form  to  M- ary  symbol  decisions  or,  especially,  to  ML 
decoding.) 
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5.  Modulation  Formats 


In  this  section,  we  will  apply  the  general  formulas  from  the  previous  sections  to  specific  modulations.  For 
the  most  part,  we  will  assume  prior  probabilities  and  then  calculate  the  “capacity”  one  could  achieve 
given  those  assumptions.  Sometimes,  we  will  find  that,  in  certain  regions  of  the  curves,  selecting  a 
different  set  of  priors  would  have  given  better  efficiency  or  throughput.  In  the  true  spirit  of  maximizing 
capacity,  one  would  choose  the  optimizing  prior  at  each  signal  level  as  long  as  the  selection  did  not 
violate  some  other  system  constraint. 

5.1  Generalized  On-Off  Keying  (OOK) 

OOK:  Preliminaries 

During  a  symbol  period,  Ts,  the  transmitter  sends  either  nothing  (corresponding  to  a  0  and  with 
probability  po)  or  a  pulse  (corresponding  to  a  1,  with  probability  p\)  such  that  an  average  of  Np  photons  is 
detected  by  the  receiver  during  the  pulse.  Note  that  this  is  a  binary >  modulation  even  when  p\  *  p 0.  The 
average  received  number  of  photons  per  symbol  is 


(73) 


Ns  —  p\  Np . 


We  thus  see  that  the  peak-to-average  ratio  of  photon  flux  (and  power)  is  \lp\.  Traditional  OOK  uses  the 
values p\  =  po  =  1/2  with  a  peak-to-average  ratio  of  2.  We  call  the  case  p,  1/2  generalized  OOK. 

We  note  that  the  pulse  can  actually  be  of  arbitrary  shape  although,  in  this  analysis,  we  assume  that  the 
pulses  do  not  overlap.  (Note  that  the  instantaneous  peak-to-average  ratio  would  change  for  non-fiat 
pulses.)  As  discussed  in  Section  2.3,  this  shape  must  be  known  for  optimum  performance  in  receiver 
types  that  need  to  separate  the  signal  from  additive  noises. 

OOK:  Noiseless  Photon-Counting  Hard  Decisions 

Whatever  the  (classical)  measurement  type,  a  hard  decision  on  whether  a  pulse  is  present  or  not  requires  a 
comparison  threshold.  The  receiver  declares  “1”  if  the  measurement  is  greater  than  this  threshold  and  “0” 
if  it  is  less.  The  threshold  must  be  chosen  to  optimize  the  desired  criterion. 

A  photon  counter  in  a  noiseless  system  has  the  easiest  task  of  all  the  OOK  hard-decision  receiver  types.  If 
it  detects  any  photons  at  all,  it  correctly  declares  “1.”  If  it  detects  none,  it  declares  “0”  and  makes  an  error 
only  if  a  pulse  was  actually  transmitted.  Thus,  with  the  Poisson  model,  we  see  that  error  probabilities  are 


P(  1 10)  =  0 

/’(Oil)  =  exp  (-  NP) . 


(74) 

(75) 
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Symbol  (bit)  error  probabilities  can  be  calculated  using  Equation  (13) 


Ps  =  P\  exp  (-  NP )  =  p\  exp  (-  Ns  !p\) . 


(76) 


This  bit-error  probability  is  plotted  in  Chart  12  for  the  symmetric  case,  p\  =  1/2. 

We  can  also  insert  Equations  (74)  and  (75)  into  Equations  (27)  and  (41)  to  find  that  capacity  and  cutoff, 
for  the  generalized  case,  are 

Cc  =  H[pl(l-exp(-Ns  / /?,))]—  ptH[exp(-Ns  / /?, )]  (77) 

Rc  =  -log[l-2p,p0(l-exp(-/Vs/2pI))]  .  (78) 


These  results  are  all  plotted  in  Charts  13-16  for  a  variety  of  p,. 

We  can  note  here  that  generalized  OOK  can  achieve  higher  and  higher  efficiency  as  the  peak-to-average 
power  ratio  increases.  On  the  other  hand,  as  this  ratio  increases,  there  are  certain  higher  data  rates  that 
become  unattainable  for  fixed  pulse  widths.  For  instance,  for  p\  =  1/32,  the  system  cannot  transmit  more 
than  about  0.2  bits  per  OOK  symbol. 

OOK:  Heterodyne  and  Homodyne  Receivers — Coherent  Hard  Decisions 

Heterodyne  and  homodyne  coherent  receivers  (again,  assuming  sampled  outputs  of  matched  filters) 
produce  a  signal,  v,  in  Gaussian  noise  with  normalized  density  function  description 

p(y  1 0)  =  fc(y)  (79) 


and 


/7(v|1)=/g(.v-V2*W/>  ), 


(80) 


where  .s  =  I  for  heterodyne  and  s  =  2  for  homodyne.  (Functions  here  were  defined  in  Section  2.2.) 
With  these,  and  basing  decisions  upon  a  threshold,  T,  we  find 
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P(\\0)  =  \p(y\0)dy 

T 


=  <D  (T) 


(81) 


T 

P(0  1 1)  =  Jp( y|  \)dy 


=  1-0  (T-yl2sNP). 


(82) 


Inserting  these  into  the  formula  for  symbol  error.  Equation  (13),  we  can  optimize  over  T  to  find  that  the 
optimum  T  satisfies 


PtP(T\\)  =  p0p(T\0)- 


(83) 


For  this  coherent  case,  the  result  is 


sNP  +  In  — 
Pi 

y[2  sNp 


(84) 


In  the  most-used  symmetric  case  of  p{)  =  p\,  we  see  that 

T  =  yfsN  p  /  2 , 


(85) 


which,  when  inserted  into  Equation  ( 13),  produces 

ph  =  O {-sjsNp  /2) 

=  0(7^). 


(86) 


Notice  that  optimum  performance  is  achieved  only  if  the  receiver  selects  the  threshold  optimally,  an 
action  which  requires  knowledge  of  the  signal  power  at  the  receiver.  A  system  that  deduces  this  level  in 
real  time  is  known  as  an  Automatic  Gain  Control  (AGC). 

These  heterodyne  and  homodyne  probabilities  for  p0  =  p\  are  plotted  in  Charts  17  and  22.  We  can  also  use 
Equations  (81),  (82),  and  (85)  in  Equations  (27)  and  (41)  to  calculate  channel  capacity  and  computation 
cutoff  rate  for  hard  decisions  with  general  p,.  Unfortunately,  the  values  of  the  threshold,  7,  that  optimize 
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capacity  and  cutoff  are  different  from  that  derived  in  Equation  (84).  These  values  are  best  calculated 
numerically.  All  these  results  are  plotted  in  Charts  18-21  and  Charts  23-26. 


OOK:  Preamplified  and  Heterodyne  Receivers — Noncoherent  Hard  Decisions 

In  the  preamplified  and  heterodyne  noncoherent  receivers,  we  use,  instead  of  the  Gaussian  density,  the 
Rician  density  (defined  in  Equation  (7): 


p(y\  0  =  /*,<■  (V2yvs,  v) 

/>(y  1 0)  =  /*„.((),  y). 


Using  a  threshold,  7,  as  in  Equations  (81)  and  (82),  we  find,  using  Equation  (10) 


P(l|0)=  0(0,7-) 

(89) 

P(0|1)=1-C9(V^>)  • 

(90) 

Inserting  these  into  the  formula  for  symbol  error.  Equation  (13),  and  optimizing  over  7,  we  find  that  7  is 
defined  by  (for  the  symmetric  case  of  p()  =  p\) 

exp {Np)  =  I0(T-j2N~p ) 

or 

T=  J — /  o~'  (ex  p(Np  », 

which  must  be  solved  numerically  although  it  has  a  good  approximate  solution  as 

T  =  y[2  +  N^J2  •  (93) 


(91) 


(92) 


Inserting  this  hack  into  Equation  (13),  we  can  calculate  the  error  probability,  which  is  shown  in  Chart  27. 
We  can  also  use  Equations  (89),  (90),  and  (92)  in  Equations  (27)  and  (41 )  to  calculate  channel  capacity  and 
computation  cutoff  rate  for  these  hard  decisions.  As  in  the  coherent  hard  decision  section,  the  thresholds 
must  be  numerically  reoptimized  for  capacity  and  cutoff.  The  results  are  plotted  in  Charts  28-3 1 . 
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OOK:  Quantum  Receiver — Hard  Decisions 

The  two-dimensional  optimization  problem  described  by  Equations  (60)  and  (61)  can  be  solved  using 
polar  coordinates  with  solutions  given  as  trigonometric  equations.  These  will  not  be  presented  here. 
However,  the  result  can  be  shown  to  provide  performance  described  by  [12]: 


P(\ |0)  =  1/2^1- 


V'  -4/>iA)  exp  (-Nr) 


[1-2/?,  exp(-A^)] 


(94) 


P(0|1)  =  1/2]  1  — 


(95) 


a/1-4Pi  Po  exp(-A^p) 


[1  -  2/?()  e\p(-N,, )]  >  , 


from  which  we  can  calculate 


Ps  =  ]/2(l-7l  -4p,p0  exp {-Np)') 


(96) 


These  have  been  plotted  in  Chart  32  for  p{)  =  p\  =  1/2  .  We  can  also  use  Equations  (94)  and  (95)  in 
Equations  (27)  and  (41)  to  calculate  semiclassical  channel  capacity  and  semiclassical  computation  cutoff 
rate  for  these  hard  decisions.  These  results  are  shown  in  Charts  33-36. 

OOK:  Photon-Counting  Soft  Decisions 

In  the  noiseless  case,  we  can  see  that  there  is  no  more  information  if  we  receive  2  or  3  or  10  photons  than 
if  we  receive  1 .  Thus,  the  soft  decision  capacity  and  cutoff  are  the  same  as  for  hard  decisions.  This  would 
not  be  the  case,  however,  if  there  were  noise  mixed  in  with  the  signal. 

OOK:  Heterodyne  and  Homodyne  Receivers — Coherent  Soft  Decisions 

To  calculate  capacity  using  coherent  soft  decisions,  we  insert  Equations  (79)  and  (80)  into  Equation  (36). 
We  find 


cc  = 


\_ 

2 


\og(2ne)~  j p{y)\og  p{y)dy 


(97) 


with 


P(  y )  =  Pofc  ( >’) +  P\  fa  ( >’ "  V 2'vyV p  )  ’ 


(98) 
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which  requires  numerical  integration.  (Remember,  s  =  1  for  heterodyne  and  .s  =  2  for  homodyne.) 
Computation  cutoff  is  similarly  found  to  be 


Rc  =  —log  [l—  2  poPl  +2p0pt(B(fG(y-^2sNP),fG(y)) 


(99) 


where  the  Bhattacharya  distance  for  Gaussian  densities  has  the  closed  form  solution 

«(/(,  (>’  -  psNp  ),  f G  (>’))  =  exp(-sN ,,  / 4)  . 


(100) 


These  results  are  plotted  in  Charts  37^)4. 

OOK:  Prcampliflcd  and  Heterodyne  Receivers — Noncoherent  Soft  Decisions 

To  calculate  capacity  using  noncoherent  soft  decisions,  we  insert  Equations  (87)  and  (88)  into  Equation  (35). 
We  find 


Cc  =  -Po  J/«„  (O.T)!og 


f  Ric  v) 


dy 


-  Pl  J )  log 


P\  +  Po 


/  Rici'J'ZN  p,y) 


dy 


Computation  cutoff  is  similarly  found  to  be 


Rc  =  -log[l-2p0p,  +  2p0pt «(/««.(>/ 2NP  ,  v),  fRu  (°,  v))], 

where  the  Bhattacharya  distance 


(101) 


(102) 


y  00  ^ 

y)’fnc (0.  v))  =  exp(-iVp  / 2)  J*  y exp(— y2  /2)yj I0(yj2N,,  )dy 


(103) 


does  not  have  a  closed-form  solution  and  so  must  be  calculated  numerically.  Capacity  and  cutoff  results 
are  plotted  in  Charts  45-48. 
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OOK:  Quantum  Receiver — Maximum-Likelihood  Decoding 

For  quantum  decoding,  we  need  to  calculate  the  eigenvalues  of  the  density  operator  of  Equation  (62).  For 
coded  OOK,  it  can  be  shown  that  [12] 


A,  =  1/2 


l±-y/l  -4p,  /?0(1  -exp(-/Vr)) 


(104) 


and  thus 

CQ  =  H(A0)=H(A])  (105) 

and 


Rq  =  -\og(Al  +A; ) 

=  -log[l-2/?lp0(l-exp(-A(p  ))]  • 


(106) 


These  results  are  plotted  in  Charts  49-52. 

5.2  Orthogonal  Modulations — PPM  and  FSK 

Orthogonal:  Preliminaries 

An  orthogonal  modulation  is  one  such  that  the  integral  of  the  product  of  any  two  different  waveforms  in 
the  set  is  zero.  Orthogonal  modulations  can  be  either  binary  or  M-ary  for  integer  M.  The  orthogonal 
dimensions  can  be  time,  frequency,  polarization,  or  other  features.  The  most  popular  time-orthogonal 
waveform  for  lasercom  is  the  pulse  position  modulation  (PPM)  set.  Here,  time  is  divided  up  into  M 
consecutive,  non-overlapping  timeslots.  Exactly  one  of  the  M  carries  a  pulse,  while  the  others  are  empty. 
Usually,  M  is  selected  to  be  a  power  of  2,  M  =  2K.  This  way,  an  M-ary  channel  symbol  represents  exactly 
K  channel  bits. 

Notice  that  orthogonal  modulations  have  the  symmetry  property  discussed  earlier.  Thus,  we  will  assume 
equiprobable  prior  probabilities,  p ,  =  I/M. 

We  note  here  that  PPM  is  very  much  like  OOK  with  p\  =  I/M.  That  is,  they  both  have  the  same  peak-to- 
average  ratio  and  the  same  probability  of  sending  a  pulse.  The  difference  is  that  PPM  requires  that  there 
be  exactly  one  pulse  every  M.  where  OOK  has  no  such  restriction. 

Frequency  shift  keying  (FSK)  has  a  single  pulse  per  symbol,  but  it  is  a  pulse  of  one  of  M  wavelengths,  or 
frequencies.  These  are  selected  so  that  the  time  integral  of  the  product  of  two  of  the  sinusoids,  over  the 
period  of  the  pulse,  is  zero.  This  property  is  true  when,  for  instance,  the  tones  are  separated  by  a 
frequency  \/TP,  where  TP  is  the  time  duration  of  the  tone  pulse. 
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We  can  also  send  a  pulse  of  one  polarization  or  the  orthogonal  polarization  (e.g.,  left-hand  circular  and 
right-hand  circular).  Although  the  modulation  called  polarization  shift  keying  (POLSK)  is  more  general, 
its  simple,  binary  orthogonal  form  is  in  the  class  of  orthogonal  modulations. 

Finally,  we  will  see  later  that  binary  DPSK  has  the  orthogonality  property,  too,  although  over  two  bits  at  a  time. 

The  reason  we  group  together  these  disparate  waveforms  is  that  they  will  be  seen  to  have  the  same 
analysis  and  performance. 

Relationships  between  pulse,  bit,  and  symbol  photon  counts  for  uncoded  orthogonal  signals  are  thus 

NP  =  NS  (107) 

=  Nh  log(A/).  (108) 

Hard  decisions  for  orthogonal  signals  are  different  from  OOK  in  that  no  threshold  is  required.  The 
receiver  first  losslessly  separates  the  M  orthogonal  dimensions — with  different  time  slots  for  PPM, 
different  band  pass  filters  for  FSK,  and  so  on.  Measurements  are  then  made  on  the  M  components.  The 
receiver  selects  the  one  that  maximizes  the  likelihood  of  choosing  the  correct  signal,  usually  the  largest 
measurement  value. 

We  will  plot  capacity  and  cutoff  as  source  bits  per  channel  usage  (i.e.,  PPM  slot)  even  though  the 
formulas  given  are  per  channel  symbol.  These  are  then  plotted  versus  average  photons  per  usage  as  well. 
These  choices  simplify  the  comparisons  with  OOK. 

When  we  plot  efficiencies  achievable  at  channel  capacity  and  cutoff  rates,  we  would  like  to  plot  them 
versus  the  total  bandwidth  expansion  that  the  coding  plus  the  modulation  requires.  The  coding  expansion 
we  know  to  be  1/C.  However,  orthogonal  modulation  dimensions  require  extra  bandwidth.  M- ary  uncoded 
PPM  time  slots,  for  M  =  2K,  require  M  times  the  bandwidth  of  a  single  slot  (channel  usage)  to  deliver  K 
bits.  Thus,  the  PPM  modulation  expansion  is  A//log2(A/).  This  value  is  then  multiplied  by  the  expansion 
due  to  the  coding.  For  instance,  uncoded  binary  orthogonal  requires  2  times  the  bandwidth  of  a  single 
slot.  Notice  that  uncoded  4-ary  also  requires  this  expansion. 

M-ary  FSK  sends  only  one  pulse  per  symbol,  but  the  M  optical  frequencies  must  be  separated  by 
multiples  of  1  /Ts.  For  M- ary,  we  have  M  tones  to  deliver  K  bits,  with  a  similar  bandwidth  expansion  to 
PPM.  (Actually,  the  expansion  is  more  like  A/-1  because  that  is  how  many  frequency  spaces  there  are 
between  M  tones.  Nevertheless,  using  the  A//log2(M)  expansion  formula  for  FSK  is  indicative  of  the 
trades  involved.  For  exact  measurements  of  the  bandwidth  expansion,  one  would  need  to  know  the  exact 
spectrum  of  the  transmitted  pulses.)  (If  the  system  designer  is  more  interested  in  the  required  speed  of 
electronics  rather  than  in  signal  bandwidth,  he  might  find  it  useful  to  normalize  all  these  charts  in  some 
other  way.) 
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Orthogonal:  Photon-Counting  Hard  Decisions 

A  photon  counter  is  used  to  make  the  M  measurements  on  the  losslessly  separated  components  of  the 
noiseless  signal.  If  any  detector  sees  a  photon,  then  the  signal  is  correctly  announced,  and  no  error  is 
made.  An  error  is  made  only  if  no  photon  is  counted.  This  is  known  as  an  erasure  because  the  receiver 
knows  it  doesn’t  know  the  answer.  Thus,  there  are  really  M  +  I  possible  outputs — M  correct  ones,  plus  an 
erasure.  Sometimes,  the  information  that  an  erasure  has  taken  place  is  useful  to  the  subsequent  signal 
processing,  and  sometimes  it  is  not.  If  it  is  not,  the  receiver  can  flip  an  M-sided  coin  to  make  a  guess. 

The  probability  of  getting  an  erasure  is  seen  to  be 

P{erasure  \  k)  =  txp(-N  p )  ( 109) 

and  we  then  have 

P(k  \k)  =  \-exp(-N p) .  (110) 

Thus,  the  probability  of  a  symbol  error,  really  an  erasure,  is 

Ps  =  exp(-Np) ,  (III) 


which  is  shown  in  Chart  53. 

If  the  decoder  does  not  allow  erasures,  then  the  probability  of  a  symbol  error,  after  the  coin  flip,  is 


ps  = 


M~\  ,  */  s 

— —  exp(-NP) 

M 


(112) 


since  it  will  map  an  erasure  into  a  correct  symbol  only  one  time  out  of  M.  Using  Equation  ( 19)  we  see  that 
such  symbol  errors  map  into  bit  errors  as 


Ph  =  ^exp(-Nh  log  A/)' 


(113) 


shown  in  Chart  54. 

We  can  calculate  capacity  and  cutoff,  assuming  the  decoder  knows  how  to  make  use  of  erasures.  The 
results  are 


Cc  =  ( 1  -  exp(-Np ))  log(M ) 

(114) 

Rc  =  log  M  -  log[l  +  (A/  -  l)exp (~NP )] . 

(115) 

Results  are  plotted  in  Charts  55-58. 
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Orthogonal:  Heterodyne  and  Homodyne  Receivers — Coherent  Hard  Decisions 

Similar  to  OOK,  we  need  to  calculate  the  probability  density  function  of  outputs  with  and  without  signals 
in  them.  Let  us  call  them  correct  and  incorrect  outputs,  with  output  densities  denoted  pc(y)>  and  p/(y).  A 
symbol  will  be  incorrectly  announced  if  any  of  the  incorrect  outputs  is  larger  than  the  correct  output. 
Equivalently,  an  error  will  be  made  if  it  is  not  true  that  the  correct  output  is  the  largest.  This  hard-decision 
symbol  error  probability  is  thus 


Ps  =  1-  JVc(y) 


j \p,(x)dx 


dy 


(116) 


For  our  heterodyne  and  homodyne  receivers,  we  have  seen  in  the  previous  section  that 


Pi(y)=fa(y) 


(117) 


and 


Pc(y)  -  /c(y-V 2sNp  )  , 


where  5  =  1  for  heterodyne  and  s  =  2  for  homodyne.  In  general,  the  integrals  in  Equation  (116)  need  to  be 
calculated  numerically,  although  for  M  =  2,  they  can  be  shown  to  simplify  to 

<*>  (V^O-  (N9) 

The  results  are  shown  in  Charts  59,  60,  65,  and  66.  Capacity  and  cutoff  can  be  calculated  by  using  this 
formula  for  Ps  and  observing  symmetries: 


Cc  =  IogM  -H(Ps)-Ps\og(M  -1) 


(120) 


Rc  =  log  M- 


2log[-v/l  —  Ps  +^(M-\)PS 


(121) 


These  results  are  plotted  in  Charts  61-64  and  Charts  67-70. 

Orthogonal:  Preamplified  and  Heterodyne  Receivers — Noncoherent  Hard  Decisions 

For  noncoherent  receivers,  we  saw  in  the  last  section  that 

Pi(y)  =  fmc(0,y) 

pc(y)=  /r„(V2a^,.v). 


(122) 

(123) 
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We  use  these  functions  in  Equation  ( 1 16)  to  find  the  symbol  error  probability 


Ps  =  1  -  jexp[—  (jc+  sNs )]  I((j4sNs  jc)  [  1  -  exp(-*)]  u  '  dx . 


(124) 


o 


The  term  in  the  square  brackets  can  be  further  expanded  and  the  M  terms  integrated  separately,  to  get  the 
final  result 


(125) 


This  formula  is  good  for  small  and  medium-sized  M.  The  formula  with  large  M ,  however,  ends  up  with 
numerical  problems,  and  then  Equation  (124)  needs  to  be  integrated  numerically. 

These  results  are  shown  in  Charts  71  and  72.  The  formulas  for  capacity  and  cutoff  are  the  same  as 
Equations  ( 1 20)  and  (121),  with  these  new  values  of  Ps  inserted.  These  results  are  shown  in  Charts  73-76. 

Orthogonal:  Quantum  Receiver — Hard  Decisions 

To  calculate  symbol  error  probabilities  using  Equation  (67),  it  can  be  shown  that  the  eigenvalues  of 
Equations  (64)  and  (65)  are  given  by 


v,  =  \+e  Ns  ^exp(-2/rt7  /  M ) 


M 


(126) 


which  simplifies  to 


(127) 


and 


(128) 


We  can  insert  these  values  into  the  formulas  for  symbol  error  probability  to  deduce 


29) 
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and  also 


(130) 


Together,  these  give  the  result 


2 


(131) 


The  symbol  and  bit-error  rates  are  shown  in  Charts  77  and  78.  The  formulas  for  semiclassical  capacity 
and  semiclassical  cutoff  are  the  same  as  Equations  (120)  and  (121),  with  this  new  value  of  Ps  inserted.  It 
can  be  shown  that  the  formula  for  cutoff  reduces  to 


Rc  =  log  M  —  log  v0  . 


(132) 


These  results  are  shown  in  Charts  79-82. 

Orthogonal:  Photon-Counting  Soft  Decisions 

As  in  the  case  of  OOK,  we  observe  that,  in  the  noiseless  case,  detecting  one  signal  photon  has  the  same 
information  as  detecting  more  than  one  signal  photon.  Thus,  soft-decision  capacity  and  cutoff  are  the 
same  as  for  hard  decisions,  with  the  understanding  that  erasures  are  kind  of  a  special  case. 

Orthogonal:  Heterodyne  and  Homodyne  Receivers — Coherent  Soft  Decisions 

To  calculate  soft-decision  performance,  we  need  to  find  the  probability  density  function  of  the  group 
(vector)  of  M  measurements.  We  can  see  that  the  y  variable  of  our  output  is  a  vector  made  up  of  { v(),  y,,  . 
...Vam  }.  Its  probability  density  function  can  be  written 


(133) 
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We  repeat  Equation  (38)  here  with  explicit  notation  of  the  vectors: 


Cc  =  log(M  )  -  lpY\x(y  I  0)  log 


A/-I 


l+I 


Py\x(y\k) 


^/^(Tl0) 


dy 


Inserting  Equation  (133)  into  Equation  (38),  we  find 


Cc  =  log(A/ )  -  JVv|*(  >'|  0)log 


M- 1 


I  +  y  Pc(y'k)  Pi  ( To ) 
it  p,(yk)  Pc  (y o) 


d  v 


(38) 


(134) 


Using  the  density  formulas  in  Equations  (117)  and  (118),  we  can  show  that  the  term  in  the  summation 
reduces  to 


exp 


(v*  -y0)xJ2sNp 


(135) 


where  the  y  variables  are  Gaussian,  the  0lh  variable  having  the  signal  as  mean.  [See  Equation  ( I  I8).| 

For  M  larger  than  2,  we  resort  to  Monte  Carlo  simulations  to  calculate  this  integral.  The  results  are  shown 
in  Charts  83  and  84,  and  in  Charts  87  and  88. 

Cutoff  rate  is  simpler.  Repeating  Equation  (48)  with  explicit  notation  of  the  vectors. 


R(  =  log (M)-  log 


2js[p(y|0),p(y  \k)\ 


_  * 

observing  the  symmetries,  and  using  Equation  ( 133),  we  can  show 


R,  =  1  og  (M  )-l  og  (l + (M - 1 W  [  pc.  Pi  ] ) . 


(136) 


(137) 


For  this  coherent  case,  we  find 

(B\pc,p,\  =  exp(-.viV,>  / 4) . 

These  results  are  shown  in  Charts  85  and  86,  and  in  Charts  89  and  90. 


(138) 
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Orthogonal:  Preamplified  and  Heterodyne  Receivers — Noncoherent  Soft  Decisions 

For  noncoherent  soft-decision  capacity,  we  follow  the  development  of  the  coherent  case,  but  use  the 
density  functions  of  Equations  (122)  and  (123).  The  summation  term  in  Equation  (134)  then  becomes  [5] 

Kiy.Wp)  _  <l39: 

/o(VoV2^  ) 


and  once  again,  we  use  a  Monte  Carlo  simulation  to  calculate  the  integral.  The  result  is  shown  in  Charts 
91  and  92. 

For  cutoff,  we  use  Equation  ( 137),  where,  for  this  noncoherent  case,  we  can  calculate 

r  i  7  / - p=-  (140 

®[  Pc,  pi  ]=  exp(-A^  /2)  y*exp(-r~ /2)y] I0(r*j2Np  )dr 

o 

[which  was  shown  in  Equation  ( 103)].  The  cutoff  result  is  shown  in  Charts  93  and  94. 

Orthogonal:  Quantum  Receiver — Maximum-Likelihood  Decoding 

For  quantum  decoding,  we  again  use  the  eigenvalues  of  the  density  operator,  which  were  shown  in 
Equations  ( 1 27)  and  (128).  We  can  insert  these  values  into  the  formulas  for  capacity  and  cutoff  to  find 

Cq  =  logAf-T  K  log  v0  +  ( M  -  1  )vt  log  vt  ]  ( 1 4 1 

M 


and 


Rq  =  log(M  ) -  log(l  +(M  - \)e~2Ns )  .  (142) 

We  can  see  that  this  is  exactly  3  dB  better  than  the  semiclassical  cutoff  rate  shown  in  Equation  ( 1 32). 

These  results  are  all  plotted  in  Charts  93-98. 
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5.3  M- ary  Phase  Shift  Keying — MPSK 
PSK:  Preliminaries 

Phase  shift  keying  uses  a  single  pulse  per  symbol,  but  the  pulse’s  optical  phase  is  one  of  M  phases, 
equally  spaced  around  the  unit  circle.  Thus,  binary  PSK  uses  0,  ft;  ternary  PSK  uses  0,  ±2ft/3;  4-ary  PSK 
(known  as  QPSK)  uses  0,  ft,  ±ft/2;  and  so  forth.  We  see  that  Ps  =  Pp  ■  We  also  note  that  the  pulses  need 
not  be  constant  amplitude.  In  fact,  any  pulse  shape  can  be  used  to  send  PSK  information,  with  peak-to- 
average  properties,  as  well  as  matched  filter  designs,  having  straightforward  extensions. 

We  can  note  that  the  MPSK  signal  is  created  in  two  dimensions  only — in-phase  and  quadrature  of  a  single 
carrier  in  a  single  timeslot.  (2PSK  uses  only  one  dimension.)  We  saw  that  orthogonal  modulations  use  M 
dimensions.  Since  only  two  dimensions  are  used  to  send  more  than  two  bits,  MPSK  for  M  >  2  is 
particularly  bandwidth  efficient.  Thus,  it  is  usually  employed  to  achieve  high  bandwidth  efficiency  rather 
than  the  highest  power  efficiency,  although  at  M  between  3  and  8,  one  can  have  some  of  both. 

The  MPSK  receiver  needs  to  deduce  which  phase  was  sent  in  each  pulse.  As  in  the  previous  modulations, 
the  best  performance  is  usually  achieved  when  the  carrier’s  phase  is  known  or  deduced  at  the  receiver 
(i.e.,  it  is  coherent  or  quantum  optimum.)  However,  noncoherent  techniques  exist  for  phase  modulations, 
too.  By  measuring  the  difference  in  phase  between  subsequent  pulses,  we  can  deduce  properties  ot  the 
phase  without  having  an  absolute  coherent  reference.  These  differential  phase  shift  keying  (DPSK) 
techniques  will  be  discussed  in  the  photon-counting  and  noncoherent  sections. 

There  is  another  interesting  facet  about  PSK  communication  that  makes  its  analysis  and  receivers 
different  from  OOK  and  orthogonal.  Since  the  possible  signals  are  not  orthogonal,  it  will  be  seen  that 
some  errors  are  more  likely  than  others.  That  is,  even  though  MPSK  has  the  symmetry  property  shown  in 
Equation  (16),  when  a  particular  phase  is  sent,  the  receiver  is  more  likely  to  make  the  error  of  choosing  a 
nearby  phase  than  a  farther  away  phase.  Soft-decision  receivers,  of  course,  implicitly  make  intelligent  use 
of  this  fact  by  measuring  the  exact  received  phase  and  then  requiring  the  decoder  to  make  the  best 
maximum-likelihood  use  of  that  information. 

Hard-decision  receivers,  for  orthogonal  signals,  are  just  as  likely  to  choose  ANY  incorrect  signal  as  any 
other.  Thus,  if  a  hard  decision  is  made,  the  decoder  for  such  a  system  performs  its  maximum-likelihood 
calculation  with  no  information  other  than  the  symbol  decision.  For  MPSK,  however,  it  a  hard  decision  is 
made,  we  see  that  we  have  two  choices.  First,  the  decoder  can  always  be  built  like  the  one  for  orthogonal 
signaling — it  just  takes  the  symbol  or  bit  decisions  and  makes  algebraic  or  other  simple  decisions.  The 
second  option  is  to  take  the  hard  symbol  decisions  and  then  require  that  the  decoder  take  into  account  the 
structure  of  the  signal  constellation.  That  is,  it  makes  a  real  maximum-likelihood  decision  on  what  it 
knows  about  which  symbols,  if  there  are  errors,  are  the  most  likely  to  have  been  substituted.  (In  reality, 
this  extra  information  requires  the  knowledge  of  the  probabilities  of  symbols  being  in  error,  and  this  is  a 
function  of  the  signal  and  noise  levels.  However,  knowledge  of  these  values  does  not  require  symbol-by¬ 
symbol  measurements  and  can  be  gotten  by  much  slower  and  simpler  AGC-type  circuits.  It  is  also  likely 
that  a  very  coarse  approximation  can  be  made  for  these  two  values  with  very  little  loss  in  performance.) 
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One  really  needs  to  weigh  the  costs  and  benefits  of  building  each  of  these  structures.  The  main  reason  one 
builds  a  simple  hard-decision  receiver  is  that  the  technology  for  making  voltage  or  other  soft-decision 
measurements  at  that  symbol  rate  is  either  unavailable  or  perhaps  just  too  expensive.  However,  even  to 
make  a  so-called  M- ary  hard  decision,  one  still  needs  to  build  a  circuit  to  compare  the  M  voltages  and 
choose  the  largest.  This  may  turn  out  to  be  almost  as  difficult  as  making  the  soft  decisions  in  the  first 
place. 

A  decoder,  too,  is  selected  not  only  to  be  paired  with  the  selected  receiver,  but  also  because  of  the 
availability  of  its  own  technology.  Maximum-likelihood  decoders  will  always  require  more  calculations 
per  bit  than  algebraic  or  other  simpler  decoders. 

Selecting  a  receiver  architecture  thus  requires  a  careful  trade  of  all  the  available  technologies.  For 
completeness,  then,  in  this  PSK  section,  we  will  present,  in  addition  to  the  soft-decision  receivers,  both 
kinds  of  hard-decision  systems:  the  true  hard  decision  receiver,  in  which  the  decoder  is  not  required  to 
make  use  of  signal  constellation  knowledge,  and  the  maximum-likelihood  decoder  based  on  the  hard 
decisions,  but  which  DOES  make  use  of  the  signal  structure.  We  will  see  that  this  latter  option  often  gives 
quite  good  performance  but,  as  we  have  stated,  may  be  an  impractical  option. 

PSK:  Photon-Counting  Hard  Decisions 

Photon  counters  cannot  discern  the  phase  of  the  carrier,  so  these  detectors  are  not  applicable  to  direct 
measurements  of  phase.  However,  there  are  at  least  two  ways  to  turn  2PSK  into  pulses  that  can  be 
detected  using  photon  counters — differentially  and  coherently.  In  Figure  8  we  show  an  architecture  for 
differential  photon  counting  of  2PSK. 


Figure  8.  Photon-counting  DPSK  receiver. 


This  receiver  splits  the  signal,  sends  the  signals  through  two  arms  whose  lengths  differ  by  one  symbol 
time,  and  then  interferes  the  two  signals.  (That  is,  it  coherently  adds  and  subtracts  the  two  signals  at  the 
output  ports.)  The  output  is  seen  to  be  the  same  as  binary  orthogonal — one  arm  or  the  other  has  all  the 
signal.  Each  output  pulse  is  made  up  of  the  sum  of  two  pulses,  each  at  half  the  power.  Thus,  this  receiver 
has  all  the  same  efficiency  properties  as  photon-counting  binary  orthogonal.  (It  does  NOT  have  the  well- 
known  3-dB  advantage  that  heterodyne  or  preamplified  DPSK  has  over  heterodyne  or  preamplified  binary 
orthogonal — as  we  will  see  soon — because  of  the  power  split.  Also,  it  should  be  noted  that  this  system 
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has  twice  the  bandwidth  efficiency  of  orthogonal  because  this  receiver  puts  out  one  bit  after  every  pulse, 
whereas  PPM,  for  instance,  needs  two  time  slots  per  bit.) 

Notice  that  this  receiver  really  produces  the  change  in  phases  between  two  bits.  Thus,  if  we  want  the 
original  bits  at  the  receiver,  we  need  to  have  performed  pre-encoding  at  the  transmitter.  The  required 
encoder  is  a  rate-one  coder  whose  output  is  the  (mod-2)  sum  of  all  the  source  bits  transmitted  so  far. 
Then,  at  the  receiver,  the  difference  between  two  phases  turns  out  to  be  an  original  source  bit.  This 
encoder  is  called,  equivalently,  a  differential  encoder  or  an  accumulator. 

There  are  at  least  two  methods  for  performing  coherent  photon  counting.  The  first  is  shown  in  Figure  9. 


Signal 

Local 
Oscillator 
Matched  in 
Phase  and 
Amplitude 

Figure  9.  Photon-counting  2PSK  receiver — balanced  LO  version. 


This  local  oscillator  is  very  special:  it  not  only  has  exactly  the  same  optical  phase  as  the  signal,  as  in  a 
homodyne  LO,  but  also  has  exactly  the  same  amplitude.  Thus,  we  see  that  all  the  power  from  both  arms 
goes  to  either  one  or  the  other  output.  This  is  another  binary  orthogonal  system  [see  Equation  (110)]  but 
this  time  with  3  dB  more  power  than  the  DPSK  version. 

An  even  better  coherent-photon-counting  method  was  invented  by  Kennedy  [12]  as  an  approximation  to 
quantum  hard  decision.  On  a  /(:  I  splitter,  with  K  approaching  infinity,  we  place  the  signal  and  a  very 
large  LO.  The  LO  is  such  that,  at  the  output  where  most  of  the  signal  emerges,  the  LO  has  the  same 
amplitude  as  the  signal.  Thus,  this  output  is  either  zero  or  twice  the  input  amplitude.  (The  waste  port  has 
almost  no  signal  and  a  very  large  amount  of  LO.)  It  therefore  has  OOK  performance  [Equation  (76)]  but 
with  4  times  the  energy  in  an  on  pulse,  which  is  3  dB  better  than  photon-counting  OOK. 

An  interesting  method  of  using  photon  counting  for  4PSK  was  proposed  by  Bondurant  [4],  It  uses  a 
configuration  like  the  Kennedy  configuration.  The  LO  starts  out  with  the  phase  of  any  of  the  four  possible 
signals.  If  the  selected  phase  turns  out  to  be  the  negative  of  the  transmitted  phase,  then  that  signal  is 
nulled  out  and  no  photons  will  be  counted  at  the  output. 
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If  the  selected  phase  is  one  of  the  three  other  phases,  then  a  signal  appears  at  the  output.  This  may  or  may 
not  lead  to  the  counting  of  a  photon.  If  a  photon  is  counted,  then  we  know  that  the  phase  of  the  LO  is  not 
the  negative  of  that  of  the  true  signal,  and  we  immediately  change  it  to  be  one  of  the  other  possible 
phases.  The  idea  is  that  we  keep  changing  phases  as  we  learn  which  signal  it  is  NOT.  Wherever  we  stop, 
either  because  we  are  right  or  because  we  ran  out  of  time,  we  announce  that  phase  as  the  negative  of  the 
actual  transmitted  signal. 

Bondurant  proposes  two  different  rules  for  deciding  on  the  subsequent  phases.  They  have  only  slightly 
different  performance,  so  we  will  analyze  only  his  so-called  Type  I  receiver. 

He  shows,  assuming  Gray  coding,  that  the  probabilities  of  counting  0  or  I  or  2  or  3  photons  in  succession, 
assuming  the  various  true  signals,  are  given  by  the  formulas  in  Table  2.  Note  that  the  columns  must  sum 
to  I. 


Table  2 

Event  Probabilities  for  the  Bondurant  Type  I  Photon-Counting  Receiver  for  QPSK 
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Here,  we  assume  that  the  user  starts  with  the  nulling  signal  for  Soo,  followed  by  50i,  SM,  and  finally  5,„, 
the  signals  walking  around  the  circle.  By  symmetry,  we  need  only  analyze  one  starting  point,  Sqo, 
especially  if  each  new  symbol  at  the  receiver  gets  a  randomly  chosen  starting  point. 

For  hard  decisions,  we  assume  that  “no  counts”  is  announced  as  Soo,  “one  count”  as  S0i.  and  so  on.  From 
this  table,  we  can  deduce  that  the  probability  of  making  a  2-bit  error  is 

[P(0|S,  ,)  +  />(  I  |S,0)] /4  =  (1/4)  exp  (-2  Ns)  (143) 

and  the  probability  of  choosing  one  bit  correct  and  one  incorrect  is 

[P(0|Sol)  +  P(0|S,o)  +  P(  I  |Sn)  +  /5(2|510)]  /  4  =  (Ns  +  1/2)  exp  (-2 Ns) .  (144) 

From  these  formulas,  we  can  calculate  symbol  and  bit-error  probabilities.  To  calculate  classical  capacity 
and  cutoff,  we  have  two  choices.  The  true  hard-decision  analysis  assumes  that  the  decoder  is  given  just 
the  symbol  decision,  with  no  information  about  the  path  it  took  to  get  there.  Thus,  the  hard-decision 
calculations  of  capacity  and  cutoff  use  only  Equations  (143)  and  (144).  We  will  examine  a  more  soft- 
decision  approach  in  a  later  section. 
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Error-rate  performances  of  all  these  (coherent)  photon-counting  PSK  systems  are  shown  in  Charts  99  and 


100.  (We  have  also  included  symbol  error  probability  for  the  Type  II  receiver,  as  calculated  by 
Bondurant.)  The  capacity  and  cutoff  results  are  shown  in  Charts  101-104.  (As  an  aside,  we  can  note  here 
that  the  cutoff  formula  for  photon-counting  orthogonal  M  =  2  and  photon-counting  OOK  p\  =  1/2  turn  out 
to  be  exactly  the  same,  as  can  be  seen  from  the  noncoherent  and  coherent/Kennedy  cutoff  curves.) 

PSK:  Heterodyne  and  Homodyne  Receivers — Coherent  Hard  Decisions 

If,  in  a  homodyne  receiver,  the  local  oscillator  has  the  phase  of  the  0  (or  the  1)  signal  of  a  binary  PSK 
input,  then  that  receiver’s  output  is  a  baseband  signal  corresponding  to  the  I’s  and  0’s,  which  can  be 
baseband  matched  filtered.  As  mentioned  in  Section  2.2,  an  optical  hybrid  can  be  used  to  make 
homodyne-like  measurements  of  QPSK.  These  are  at  baseband,  but  do  not  achieve  the  3-dB  improvement 
the  binary  homodyne  receiver  sees.  Finally,  a  heterodyne  receiver  uses  a  heterodyne  front  end,  and  then 
follows  it  with  a  classical  (RF-style)  coherent  MPSK  matched-filter  receiver. 

It  can  be  shown  that  the  hard-decision  error  probability  for  homodyne  2PSK  is 


Ph  =  ■ 


(145) 


For  coherent  heterodyne,  MPSK  has  the  symbol  probability 


(146) 


where  we  let  PSm  be  the  M-ary  symbol  error  probability,  and  where  Craig’s  function  fc  was  defined  in 
Equation  (6).  It  has  the  properties  that 


(147) 


P52  =  0(V2M~)  M  =  2 


PS4  =  \-0-PS2)2  - 


(148) 


Interestingly,  the  bit-error  probability  from  Gray-coded  4PSK  is  the  same  as  the  bit  (symbol)  error 
probability  from  2PSK. 

For  the  true  hard-decision  receiver,  as  discussed  in  the  introductory  PSK  section,  we  need  to  analyze  a 
decoder  that  assumes  only  that  a  symbol  is  either  correct  or  incorrect.  Thus,  the  error  probabilities  for 
each  of  the  (M  -  1)  incorrect  symbols  is  assumed  to  be 


Pe=Psm/(M-  1) 


(149) 


and  the  probability  of  the  correct  symbol  choice  is 


51 


(150) 


p  =  |  -  p 
rc  1  Cs 

Symbol  error  probabilities  are  shown  in  Chart  105  for  heterodyne  and  Chart  1 13  for  homodyne.  Bit-error 
probabilities  are  shown  in  Charts  106  and  1 14. 

Using  these  formulas,  we  can  also  calculate  the  pure  hard-decision  coded  performance.  We  only  show 
cutoff  for  this  case  because  any  decoding  technique  that  could  have  achieved  capacity  would  have  used  a 
maximum-likelihood  algorithm,  for  which  the  added  complexity  of  signal  structure  knowledge  would 
certainly  have  been  worth  it.  We  show  cutoff  results  for  this  heterodyne  MPSK  case  in  Charts  107  and  108. 

To  calculate  capacity  and  cutoff  for  hard  decisions  with  maximum-likelihood  decoding,  we  need  to 
calculate  the  probabilities,  not  just  of  symbol  errors,  but  of  choosing  each  output  symbol  given  each 
possible  input  symbol.  By  the  equal  phase-spacing  property,  we  see  that  we  need  only  find  the 
probabilities  of  declaring  each  k  *  0,  given  that  0  was  sent.  Let  us  number  the  symbols  0,  1, .  .  . ,  M  -  I  as 
we  travel  around  the  unit  circle.  That  is,  symbol  k  is  declared  if  the  phase  is  measured  as  being  between 

(2k  -\)n  .  (2k  +  \)n 

-  and  - . 

M  M 


Using  Craig’s  probability  function  [Equation  (6)],  it  can  be  shown  that  the  probability  of  choosing  symbol  k  is 


P(k  |  0)  =  - 
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/<■  a/2v7. 
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(151) 


=  P(M -k  10) 


and,  for  M  even,  the  center  (M/2)  term  is  given  by 


2  M 


(152) 


We  first  use  these  formulas  to  calculate  the  heterodyne  bit-error  probabilities  by  weighting  each  multi-bit¬ 
error  pattern  with  its  probability.  These  probabilities  are  shown  in  Chart  106. 

Using  these  probability  formulas  in  Equations  (28)  and  (42),  we  can  also  calculate  the  capacity  and  cutoff 
curves,  whose  results  are  shown  in  Charts  109-112  for  heterodyne  and  Charts  1 15-1 18  for  binary 
homodyne. 

PSK:  Preamplified  and  Heterodyne  Receivers — Noncoherent  Hard  Decisions 

As  mentioned  briefly  above,  the  noncoherent  techniques  for  MPSK  measure  differences  in  phase  between 
subsequent  symbols.  (Although  we  will  not  pursue  them  here,  there  are  more  efficient  noncoherent  MPSK 
receivers  that  examine  more  than  two  symbols  in  a  row  [9,  16].)  A  particularly  efficient  architecture  for 
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binary  systems  that  has  found  much  utility  in  wideband  systems  is  an  optical  preamplifier,  followed  by 
the  interferometer-based  demodulator  shown  in  Figure  8.  As  discussed  earlier,  such  a  DPSK  system 
requires  a  precoding  at  the  transmitter  with  a  differential  encoder. 

For  larger  M,  one  would  need  to  make  phase  measurements  in  the  two  symbols  and  then  deduce  which 
phase  transition  was  most  likely,  although  there  are  sometimes  architectural  tricks  to  simplify  this 
structure. 

The  symbol  error  probability  of  a  differential  MPSK  receiver  (sometimes  called  MDPSK  or  some  other 
variant)  can  be  shown  to  be 


n  =/,<V W.jj). 


(153) 


where  fP  is  Pawula’s  function,  defined  in  Equation  (11).  This  symbol  error  probability  is  plotted  in  Chart 
1 19  with  the  related.  Gray-coded  bit-error  probability  in  Chart  120. 


As  in  the  coherent  case,  we  first  calculate  the  true  hard-decision  cutoff  using  this  symbol  error  probability 
along  with  an  analog  of  Equations  ( 149)  and  (150).  These  results  are  shown  in  Charts  121  and  122. 


Also  following  the  coherent  case,  in  order  to  calculate  capacity  and  cutoff  for  hard  decisions  with 
maximum-likelihood  decoding,  we  need  to  find  specific  symbol  error  probabilities.  Exactly  paralleling 
the  coherent  case,  we  find  that  the  probability  of  choosing  symbol  k  when  symbol  0  is  correct  is 


m|o>=- 


I  M  J  \  M 


\<k< 


M-  1 


=  P(M-k  |0) 


(154) 


and  for  M  even,  we  have 


P(^-\0)=My[2N^, 


( M  -  I  )x 
M 


(155) 


Using  these  functions,  we  can  calculate  capacity  and  cutoff,  with  the  results  shown  in  Charts  123-126. 

We  can  note  here  briefly  that  binary  DPSK  compares  the  sequences  (1,1)  and  (1,-1),  which  are  seen  to 
be  orthogonal  two-bit  sequences.  Thus,  all  the  performance  equations  (in  heterodyne,  homodyne,  and 
preamplified  receivers)  for  binary  DPSK  are  the  same  as  binary  orthogonal  except  that  the  DPSK  decision 
uses  the  energy  from  two  sequential  symbols.  Thus,  binary  DPSK  is  exactly  3  dB  better  than  binary 
orthogonal  in  every  metric  (except  the  noncoherent  photon-counting  variant,  as  mentioned  earlier). 
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PSK:  Quantum  Receiver — Hard  Decisions 

Without  performing  the  derivation  here,  it  can  be  shown  that  the  optimized  quantum  measurement 
operators  result  in  performance  described  by  [  12] 


P(.i\k)  = 


1  M- 1  _ 

—  exp[~jx2(j-k)i/M  ] 


=  P(j~k  |0) 


where 


(156) 


v,  =  XexP[_2yVv 

k=  0 


\MJ 


IVl  J 


(157) 


Note  that  both  these  formulas  involve  discrete  Fourier  transforms.  Also  note  the  explicit  satisfaction  of 
the  symmetry  property  in  Equation  (156). 

Using  these  formulas,  we  can  calculate  the  symbol  probability  using  Equation  (67).  Results  are  shown  in 
Charts  1 27  and  128. 

We  reason  that  one  who  works  hard  enough  to  build  a  quantum  optimum  receiver  would  not  settle  for  the 
“pure”  hard-decision  decoder.  Thus,  we  will  only  calculate  the  maximum-likelihood  decoding 
performance.  Using  Equation  (156)  in  Equations  (28)  and  (42),  we  can  calculate  semiclassical  capacity 
and  semiclassical  cutoff  for  quantum  hard  decisions  (for  the  L  =  M  case,  as  discussed  earlier).  Care  must 
be  taken  with  the  complex  nature  of  the  xjk  values  defined  in  Equations  (55)  through  (59). 

These  results  are  all  shown  in  Charts  129-132. 

PSK:  Photon-Counting  Soft  Decisions 

As  in  the  cases  of  OOK  and  orthogonal,  we  observe  that,  in  the  noiseless  binary  cases,  detecting  one 
signal  photon  has  the  same  information  as  detecting  more  than  one  signal  photon.  Thus,  soft-decision 
capacity  and  cutoff  are  the  same  as  for  hard  decisions.  This  applies  to  all  the  binary  photon-counting 
methods  presented  earlier. 

The  Bondurant  photon-counting  4PSK  approach,  however,  has  a  kind  of  soft  metric  that  comes  along 
with  the  hard  decision — namely,  how  many  counts  were  counted.  We  would  know,  for  example,  that  if 
we  announced  “01”  after  one  count,  the  answer  is  likely  01,  but  might  also  be  1 1  or  10.  However,  it  could 
NOT  be  00,  since  that  hypothesis  was  rejected  by  counting  the  first  photon.  Thus,  the  decoder  can  make 
use  of  this  information. 
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The  symbol  error  probabilities  were  given  in  Table  2  earlier.  Assuming  that  each  starting  phase  is 
randomly  chosen,  we  can  calculate  classical  capacity  for  this  hard-decision-plus-ML  decoding  system. 
The  results  are  shown  in  Charts  133  and  134. 


We  should  also  note  that  Bondurant  designed  this  receiver  for  uncoded,  hard-decision  use.  Hence,  its 
performance  is  less  good  for  low  code  rates.  It  is  likely  that  a  variant  on  his  decision  mle  could  give  better 
coded  performance. 

There  is  one  more  observation  to  be  made  here.  It  has  been  shown  [2]  that  the  likelihood  function  for  any 
photon-counting  scheme  is  more  correctly  a  function  of  the  actual  times  of  the  photon  arrivals.  In  fact,  in 
his  Type  II,  Bondurant  used  these  times  to  try  to  improve  QPSK  hard-decision  performance.  Using  either 
Type  I  or  Type  II,  the  exact  likelihood  functions  for  receiving  0  or  I  or  2  photons  in  succession  can  be 
calculated,  using  the  arrival  times.  We  will  not  perform  this  calculation  here. 

PSK:  Heterodyne  and  Honiodyne  Receivers — Coherent  Soft  Decisions 

To  calculate  capacity  and  cutoff  rate  for  coherent,  soft-decision  PSK,  we  need  to  find  the  distances 
between  the  signal  points  around  the  unit  circle.  These  can  be  found  using  straightforward  geometry.  The 
results,  when  inserted  into  the  formulas  of  Equations  (38)  and  (48)  can  be  shown  to  be 


Q  =  log(M )  - 

{/gU)J/(;(  V)  log 


(158) 


/ 

M- 1 

Xexp 

-2< 

k= 0 

\ 

N  sin2 


7tk ^ 


K  sin  A 

M 


.  7Tk  Kk 

jcstn( — )-  ycos( — ) 
M  M 


dvdx 


and 


r  m -i 

Rc  =  log  (M )  -  logj  1  +  ^exp 

{  k= i 

Since  the  capacity  is  a  simple  (two-dimensional)  expected  value,  it  can  be  calculated  using  either 
numerical  integration  or  Monte  Carlo  techniques. 

These  results  are  shown  in  Charts  135-138. 

For  homodyne  2-PSK.  replace  Ns  by  2 Ns  in  these  formulas.  For  other  M,  homodyne  is  not  generally 
attempted.  (See  the  discussion  in  Section  2.2.) 

The  results  for  binary  homodyne  are  shown  in  Charts  139-142. 


«/  •  i 

-  N,  sin  ( — 
M 


(159) 
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PSK:  Preamplified  and  Heterodyne  Receivers — Noncoherent  Soft  Decisions 

To  calculate  the  performance  of  noncoherent  MPSK,  we  need  to  find  the  probability  density  of  the 
complex  components  of  the  two  adjacent  symbol  times,  but  averaged  over  the  unknown  phase  angle, 
assumed  to  be  uniformly  distributed.  This  is  kind  of  a  generalization  of  the  Rician  density  [Equation  (7)] 
and  is  given  by 


p{  y  |  A: )  =— — exp 
2/r 


■{in +IzoJi±1aJ1) 

2 


y0 


+  y,£ 


jk In  I  M 


(160) 


where  y  is  the  two-dimensional  complex  vector  of  the  two  symbols.  We  can  see  the  phase  relationship 
between  the  two  symbols  is  2nk/M,  which  implies  a  differential  encoding.  We  insert  this  formula  in 
Equation  (38)  to  find  the  capacity.  Similarly,  we  insert  it  into  Equation  (48)  to  find  the  cutoff.  These 
results  are  shown  in  Charts  143-146.  As  mentioned  above,  there  are  more  efficient  receivers  that  make 
use  of  multiple  sequential  symbols  [16].  Also  as  mentioned  earlier,  binary  (noncoherent,  2-bit  observation 
interval)  DPSK  has  exactly  the  same  form  as  binary  orthogonal,  but  with  3  dB  better  performance. 

PSK:  Quantum  Receiver — Maximum-Likelihood  Decoding 

We  can  use  the  eigenvalues  of  Equation  (157)  in  Equation  (66)  to  calculate  quantum  capacity.  We  can 
also  use  them  in  Equation  (70)  to  find  quantum  cutoff  explicitly  (using  Parseval’s  theorem  to  simplify  the 
formula) 


M- 1 

Rq  =  log(M)  —  log(l  4-  ^exp[-4 Ns  sin2 

k= I 

which  is  seen  to  be  exactly  3  dB  better  than  cutoff  using  homodyne  soft  decision.  These  results  are  shown 
in  Charts  147-150. 
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6.  Discussion  and  Comparison  of  Some  Modulation 
Formats  and  Receivers 


In  this  section,  we  will  make  explicit  comparisons  of  some  of  the  results  we  have  shown  in  order  to  see 
the  benefits  and  drawbacks  of  employing  the  different  methods. 

First,  in  Charts  151  and  152,  we  compare  the  theoretical  optimum  coded  performance  for  unconstrained 
classical  and  quantum  systems.  At  low  code  rates  (high  bandwidth,  low  signal  power)  there  is  a  3-dB 
jump  from  classical  cutoff  to  classical  capacity,  then  another  3  dB  to  quantum  cutoff.  Quantum  optimum 
capacity  is  seen  to  have  a  very  different  quality  in  that  it  continues  to  improve  as  the  bandwidth  increases. 

In  Chart  153,  we  plot  BER  for  symmetric  OOK  using  all  the  different  receiver  structures  we  have 
investigated.  We  can  see  that  in  uncoded  systems  at  high  signal  powers  where  the  BER  is  low,  there  is 
about  a  3-dB  gain  going  from  heterodyne  to  homodyne  receivers  and  another  3  dB  going  to  either  photon¬ 
counting  or  quantum  optimum.  For  hard-decision  coded  systems,  however,  where  we  would  operate  at 
high  BER,  the  efficiency  gains  are  not  so  clear-cut. 

In  Charts  154  and  155,  we  can  see  the  BER  for  orthogonal  signals  with  various  receivers,  for  M  =  2  and  4. 

In  Chart  156,  we  see  the  BER  for  all  the  receivers  for  2PSK,  including  the  several  photon-counting 
variants.  Once  again,  optimum  quantum  hard-decision  and  (coherent)  photon-counting  variants  do  the 
best  at  high  signal  powers. 

In  Chart  157,  we  compare  many  different  binary  coded  systems.  It  is  quite  clear  where  the  gains  are. 

In  Charts  158-162,  we  compare  the  performance  at  capacity  of  generalized  OOK  for  the  different 
receivers  at  various  peak-to-average  ratios.  Interestingly,  for  high  ratios,  noiseless  photon  counting  does 
just  about  as  well  as  quantum  optimum  hard  decision.  We  can  also  see  that  at  fairly  high  bandwidths  (i.e., 
low  FEC  code  rates),  this  family  essentially  achieves  optimum  hard-decision  capacity  performance. 

In  Charts  163-167,  we  compare  the  performance  at  capacity  of  orthogonal  modulation  for  various 
alphabet  sizes,  M.  Here,  it  is  interesting  to  see  that  at  high  M,  photon  counting  does  better  than  the  so- 
called  quantum  optimum  hard  decision.  This  is  because  the  hard-decision  system  had  been  constrained  to 
make  a  decision  between  the  M  possibilities,  while  our  photon-counting  receiver  was  allowed  an  M  +  Ith 
output,  the  erasure.  Although  we  will  not  derive  it  here,  an  optimum  [choose  L  >  M,  with  the  notation  of 
Equation  (57)]  symbol  measurement  would  do  better  than  this  L  =  M  measurement,  although  tor  high  M, 
it  is  doubtful  that  the  optimum  would  be  much  better  than  the  photon-counting  receiver. 

In  Charts  168-171,  we  compare  generalized  OOK  and  orthogonal  for  several  receiver  types.  We  can  see 
that  OOK  is  quite  a  bit  better  at  low  peak-to-average  ratios  (M),  but  that  the  differences  become  less  at 
high  M  and  high  bandwidth  (low  code  rate).  In  general,  we  can  see  that  these  pulsed  formats  are  very  near 
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their  theoretical  optimum  channel  capacities  for  many  M  in  heterodyne  (and  homodyne)  soft-decision  and 
quantum  optimum  systems. 

In  Charts  172  and  173,  we  compare  coded  performance  of  2-ary  and  4-ary  PSK  for  various  receivers.  We 
can  see  that  there  is  a  big  gap  between  all  the  symbol-by-symbol  hard-decision  systems  and  the  quantum 
optimum  system.  The  potential  gain  here  certainly  justifies  investigation  into  code  and  receiver 
construction  for  the  quantum  optimum. 

In  Chart  174,  we  compare  optimum  quantum  coded  performance  for  many  of  the  modulation  formats. 

For  completeness,  we  show  in  Charts  175-177  the  forms  of  the  eigenvalues  of  the  quantum  density 
operators  [see  Equations  (63)  to  (65)].  The  eigenvalues  for  PSK  have  a  quite  different  shape  from  those  of 
OOK  and  orthogonal. 

We  conclude  with  Table  3  (variants  of  which  can  be  found  in  many  references),  which  gathers  the  bit¬ 
error  rate  performance  formulas  for  all  the  binary  formats  we  have  presented,  using  each  of  the  receivers 
we  have  analyzed.  Although  it  is  of  interest,  we  have  seen  that  the  more  complex  channel  capacity 
formulas  are  really  the  ones  that  show  the  performance  we  seek. 
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Table  3 

Bit-Error  Rate  Formulas  for  Binary  Modulations 


Modulation  Format 

Receiver  Type 

Bit-Error  Probability 

Comments 

OOK 

Photon  Counting 

~exp  (~2Nh) 

Only  error  type  is  by 
missing  aT 

Heterodyne  Coherent 

■hJK) 

Requires  threshold 

Homodyne  Coherent 

0(V2  Nh) 

Requires  threshold 

Heterodyne 

Noncoherent 

^[\  +  Q(0,T)-Q(2yfN^,T)] 

Requires  threshold,  T 
[see  Equation  (92)] 

Quantum  Optimum 

■^■(1  --y/l  -  exp  (~2Nh)) 

Orthogonal 

Photon  Counting 

-•exp  (—Nh) 

Error  from  coin  flip 
when  erasure 

Heterodyne  Coherent 

Homodyne  Coherent 

o  ipNh ) 

Heterodyne 

Noncoherent 

^exp  (~Nh  12) 

Quantum  Optimum 

d-(l  - ^/l  —  exp  (-2 Nh)) 

PSK 

Photon  Counting 
Noncoherent 

—exp  (—Nh) 

Error  from  coin  flip 
when  erasure 

Photon-Counting 
Coherent  Symmetric 

~exP  (~2Nh) 

Error  from  coin  flip 
when  erasure 

Photon  Counting 
Coherent  Kennedy 

^-exp(-4  N„) 

Only  error  type  is  by 
missing  aT 

Heterodyne  Coherent 

O  (pNh ) 

Homodyne  Coherent 

(j4Nh ) 

Heterodyne 

Noncoherent 

•^-exp  (~Nh) 

Quantum  Optimum 

~(1  ->/l  -  exp  (“4A^(,)) 

7.  Conclusion 


It  is  hoped  that  this  collection  might  serve  several  functions: 

•  Act  as  a  single  place  for  finding  many  relevant  equations  for  many  communications  analyses 

•  Serve  as  a  reminder  of  all  the  design  possibilities  there  are  for  the  lasercom  system  designer 

•  Act  as  a  repository  of  basic  references,  explicit  and  implicit,  that  may  serve  as  a  jumping-off 
point  for  the  interested  reader 

However,  the  main  hope  of  the  author  is  that  these  side-by-side  presentations  and  comparisons  will 
motivate  the  reader  to  invent  efficient  systems  that  can  achieve  the  theoretical  possibilities. 
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Chart  5.  Classical  cutoff:  BSC  (hard  decision). 
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Chart  6.  Classical  cutoff:  Coherent  Gaussian. 
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Chart  1 1 .  Efficiency  at  quantum  cutoff:  Gaussian. 
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Chart  12.  Bit-error  probability:  00K  (p0=Pi  =  1/2)  photon  counting. 
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Chart  13.  Classical  capacity:  00K  photon  counting. 
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Chart  14.  Efficiency  at  classical  capacity:  OOK  photon  counting. 
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Chart  17.  Symbol  error  probability:  OOK  (p0  -p i  =  1/2)  heterodyne  or  preamplified  coherent. 
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Chart  18.  Classical  capacity:  00K  heterodyne  or  preamplified  coherent  hard  decision. 
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Chart  20.  Classical  cutoff:  00K  heterodyne  or  preamplified  coherent  hard  decision. 
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Chart  21.  Efficiency  at  classical  cutoff:  OOK  heterodyne  or  preamplified  coherent  hard  decision. 
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Chart  22.  Symbol  error  probability:  00K  (p{>  =  pi  =  1/2)  homodyne  coherent. 
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Chart  29.  Efficiency  at  classical  capacity:  00K  heterodyne  or  preamplified  noncoherent  hard  decision. 
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Chart  30.  Classical  cutoff:  00K  heterodyne  or  preamplified  noncoherent  hard  decision. 
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Chart  31.  Efficiency  at  classical  cutoff:  OOK  heterodyne  or  preamplified  noncoherent  hard  decision. 
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Chart  33.  Classical  capacity:  00K  quantum  hard  decision  (semiclassical). 
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Chart  35.  Classical  cutoff:  00K  quantum  hard  decision  (semiclassical). 
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Chart  38.  Efficiency  at  classical  capacity:  OOK  heterodyne  or  preamplified  coherent  soft  decision. 
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Chart  41.  Classical  capacity:  00K  homodyne  coherent  soft  decision. 
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Chart  45.  Classical  capacity:  00K  heterodyne  or  preamplified  noncoherent  soft  decision. 
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Chart  47.  Classical  cutoff:  OOK  heterodyne  or  preamplified  noncoherent  soft  decision. 


15.0 


Chart  48.  Efficiency  at  classical  cutoff:  OOK  heterodyne  or  preamplified  noncoherent  soft  decision. 
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Chart  54.  Bit-error  probability  (erasure  with  coin  flip):  orthogonal  photon  counting. 
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Chart  55.  Classical  capacity:  orthogonal  photon  counting. 
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Chart  57.  Classical  cutoff:  orthogonal  photon  counting. 
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Chart  58.  Efficiency  at  classical  cutoff:  orthogonal  photon  counting. 
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Chart  59.  Symbol  error  probability:  orthogonal  heterodyne  or  preamplified  coherent. 
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Chart  61.  Classical  capacity:  orthogonal  heterodyne  or  preamplified  coherent  hard  decision. 
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Chart  63.  Classical  cutoff:  orthogonal  heterodyne  or  preamplified  coherent  hard  decision. 
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Chart  64.  Efficiency  at  classical  cutoff:  orthogonal  heterodyne  or  preamplified  coherent  hard  decision. 


00+3’ L 


Ainjqeqojd  JOJjg  iig 


Chart  65.  Symbol  error  probability:  orthogonal  homodyne  coherent. 
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Chart  66.  Bit-error  probability:  orthogonal  homodyne  coherent. 


Chart  68.  Efficiency  at  classical  capacity:  orthogonal  homodyne  coherent  hard  decision. 
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Chart  69.  Classical  cutoff:  orthogonal  homodyne  coherent  hard  decision. 
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Chart  71.  Symbol  error  probability:  orthogonal  heterodyne  or  preamplified  noncoherent. 
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Chart  73.  Classical  capacity:  orthogonal  heterodyne  or  preamplified  noncoherent  hard  decision. 
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Chart  75.  Classical  cutoff:  orthogonal  heterodyne  or  preamplified  noncoherent  hard  decision. 
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Chart  76.  Efficiency  at  classical  cutoff:  orthogonal  heterodyne  or  preamplified  noncoherent  hard  decision. 
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Chart  77.  Symbol  error  probability:  orthogonal  quantum  hard  decision  (semiclassical). 
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Chart  78.  Bit-error  probability:  orthogonal  quantum  hard  decision  (semiclassical). 
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Chart  79.  Classical  capacity:  orthogonal  quantum  hard  decision  (semiclassical). 
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Chart  80.  Efficiency  at  classical  capacity:  orthogonal  quantum  hard  decision  (semiclassical). 


Chart  81.  Classical  cutoff:  orthogonal  quantum  hard  decision  (semiclassical). 
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Chart  83.  Classical  capacity:  orthogonal  heterodyne  or  preamplified  coherent  soft  decision. 
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Chart  87.  Classical  capacity:  orthogonal  homodyne  coherent  soft  decision. 
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Chart  89.  Classical  cutoff:  orthogonal  homodyne  coherent  soft  decision. 
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Chart  90.  Efficiency  at  classical  cutoff:  orthogonal  homodyne  coherent  soft  decision. 
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Chart  91.  Classical  capacity:  orthogonal  heterodyne  or  preamplified  noncoherent  soft  decision. 
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Chart  93.  Classical  cutoff:  orthogonal  heterodyne  or  preamplified  noncoherent  soft  decision. 
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Chart  98.  Efficiency  at  quantum  cutoff:  orthogonal. 
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Chart  99.  Symbol  error  probability:  MPSK  photon  counting. 
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Chart  101.  Classical  capacity:  MPSK  photon  counting  hard  decision. 
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Chart  103.  Classical  cutoff:  MPSK  photon  counting  hard  decision. 
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Chart  104.  Efficiency  at  classical  cutoff:  MPSK  photon  counting  hard  decision. 
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Chart  105.  Symbol  error  probability:  MPSK  heterodyne  or  preamplified  coherent. 
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Chart  107.  Classical  cutoff:  MPSK  heterodyne  or  preamplified  coherent  pure  hard  decision. 
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Chart  108.  Efficiency  at  classical  cutoff:  MPSK  heterodyne  or  preamplified  coherent  pure  hard  decision. 
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Chart  109.  Classical  capacity:  MPSK  heterodyne  or  preamplified  coherent  hard  decision  with  ML  decoding. 


1.E+01 


CD 

TJ 


o 

z 

9 


LU 


O 

a 


l 


(/) 

Q) 

C 

c 

co 


O 

7* 

c 

o 

o 

-C 

CL 


(loqiuAs  |9uueqo/s;jq  aojnos)  wojno 


Chart  111.  Classical  cutoff:  MPSK  heterodyne  or  preamplified  coherent  hard  decision  with  ML  decoding. 
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Chart  112.  Efficiency  at  classical  cutoff:  MPSK  heterodyne  or  preamplified  coherent  hard  decision 
with  ML  decoding. 
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Chart  113.  Symbol  error  probability:  MPSK  homodyne  coherent. 
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Chart  114.  Bit-error  probability:  MPSK  homodyne  coherent. 
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Chart  115.  Classical  capacity:  MPSK  homodyne  coherent  hard  decision. 
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Chart  119.  Symbol  error  probability:  MPSK  heterodyne  or  preamplified  noncoherent. 
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Chart  122.  Efficiency  at  classical  cutoff:  MPSK  heterodyne  or  preamplified  noncoherent  hard 
decision  pure  hard  decision. 
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Chart  123.  Classical  capacity:  MPSK  heterodyne  or  preamplified  noncoherent  hard  decision 
with  ML  decoding. 


15.0 


>* 

>» 

>s 

>» 

CD 

CD 

CD 

CD 

CN 

CO 

4 

CO 

1 

i 

1 

1 

i 

i 

1 

o 

4 


o 

CO 


o 

csi 


o  o  © 

r-'  ©  0> 


o 

CO 


o  o  o 

h-  CO  Lf> 


o  o 

t-  © 


o 

o 

o 


o 

o 


<D 

CL 

(ft 

<D 

U) 
TO 

V) 

z> 

0) 

c 

c 

(0 

£ 

o 


x: 

■g 

i 

TJ 

c 

CD 

DO 


(9P)  \\9  J9d  suojoqd 


o 

4 


o 

CO 


o 

CNJ 


Chart  124.  Efficiency  at  classical  capacity:  MPSK  heterodyne  or  preamplified  noncoherent  hard  decision 
with  ML  decoding. 
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Chart  125.  Classical  cutoff:  MPSK  heterodyne  or  preamplified  noncoherent  hard  decision  with  ML  decoding. 
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Chart  126.  Efficiency  at  classical  cutoff:  MPSK  heterodyne  or  preamplified  noncoherent  hard  decision 
with  ML  decoding. 
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Chart  130.  Efficiency  at  classical  capacity:  MPSK  quantum  hard  decision  (semiclassical)  with  ML  decoding. 
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Chart  131.  Classical  cutoff:  MPSK  quantum  hard  decision  (semiclassical)  with  ML  decoding. 
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Chart  132.  Efficiency  at  classical  cutoff:  MPSK  quantum  hard  decision  (semiclassical)  with  ML  decoding. 
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Chart  133.  Classical  capacity:  4PSK  photon  counting  hard  decision  with  ML  decoding. 
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Chart  134.  Efficiency  at  classical  capacity:  4PSK  photon  counting  hard  decision  with  ML  decoding. 
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Chart  135.  Classical  capacity:  MPSK  heterodyne  or  preamplified  coherent  soft  decision. 
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Chart  137.  Classical  cutoff:  MPSK  heterodyne  or  preamplified  coherent  soft  decision. 
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Chart  139.  Classical  capacity:  MPSK  homodyne  coherent  soft  decision. 
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Chart  140.  Efficiency  at  classical  capacity:  MPSK  homodyne  coherent  soft  decision. 
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Chart  141.  Classical  cutoff:  MPSK  homodyne  coherent  soft  decision. 
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Chart  143.  Classical  capacity:  MPSK  heterodyne  or  preamplified  noncoherent  soft  decision. 
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Chart  144.  Efficiency  at  classical  capacity:  MPSK  heterodyne  or  preamplified  noncoherent  soft  decision. 
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Chart  145.  Classical  cutoff:  MPSK  heterodyne  or  preamplified  noncoherent  soft  decision. 
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Chart  147.  Quantum  capacity:  MPSK. 
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Chart  149.  Quantum  cutoff:  MPSK. 
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Chart  152.  Comparison:  Efficiencies  at  classical  Gaussian  capacity,  classical  Gaussian  cutoff,  quantum 
Gaussian  capacity,  quantum  Gaussian  cutoff. 
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Chart  153.  Comparison:  BER,  00K  (p\  =  1/2) — preamplified  noncoherent,  heterodyne  or  preamplified 
coherent,  homodyne  coherent,  photon  counting,  quantum. 
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Chart  154.  Comparison:  BER,  orthogonal  (M  =  2) — preamplified  noncoherent,  heterodyne  or  preamplified 
coherent,  homodyne  coherent,  photon  counting,  quantum. 


1.E+00 


CM 

CO 

w 

ID 

h- 

CO 

0) 

O 

CN  • 

O 

o 

o 

© 

O 

o 

o 

o 

T— 

y— 

UJ 

UJ 

UJ 

UJ 

UJ 

LU 

UJ 

LU 

UJ 

UJ 

UJ 

T- 

y— 

y- 

T~ 

T— 

y- 

T— 

y— 

r-’ 

T“ 

r-’ 

CO 

■O 


CO 

V- 

5 

Q. 

</> 

c 

o 

o 

£ 

CL 


A^uqeqojd  JOJJ3  \\q 


Chart  155.  Comparison:  BER,  orthogonal  (M  =  4) — preamplified  noncoherent,  heterodyne  or.  preamplified 
coherent,  homodyne  coherent,  photon  counting,  quantum. 
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Chart  156.  Comparison:  BER,  PSK  (M  =  2) — preamplified  noncoherent  (DPSK),  heterodyne  or  preamplified 
coherent,  homodyne  coherent,  photon  counting,  quantum. 
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Chart  157.  Comparison:  Efficiency  at  capacity,  binary  modulations — OOK  (p(  =  1/2),  PPM,  and  PSK,  each  at 
preamplified  noncoherent,  heterodyne  or  preamplified  coherent,  and  quantum,  or  Gaussian  quantum. 


Chart  158.  Comparison:  Efficiency  at  capacity,  OOK,  p\  =  1/2,  soft  and  hard  decisions,  preamplified 
noncoherent,  heterodyne  or  preamplified  coherent,  homodyne  coherent,  photon  counting,  quantum,  or 
Gaussian  quantum. 
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Chart  159.  Comparison:  Efficiency  at  capacity,  OOK,  px  =  1/4,  soft  and  hard  decisions,  preamplified  noncoherent, 
heterodyne  or  preamplified  coherent,  homodyne  coherent,  photon  counting,  quantum,  or  Gaussian  quantum. 
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noncoherent,  heterodyne  or  preamplified  coherent,  homodyne  coherent,  photon  counting,  quantum,  or 
Gaussian  quantum. 
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Chart  161.  Comparison:  Efficiency  at  capacity,  OOK,  p\  =  1/64,  soft  and  hard  decisions,  preamplified 
noncoherent,  heterodyne  or  preamplified  coherent,  homodyne  coherent,  photon  counting,  quantum,  or 
Gaussian  quantum. 
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Chart  163.  Comparison:  Efficiency  at  capacity,  orthogonal,  M  =  2,  soft  and  hard  decisions,  preamplified 
noncoherent,  heterodyne  or  preamplified  coherent,  homodyne  coherent,  photon  counting,  quantum,  or 
Gaussian  quantum. 


o 

o 

o 


o 

o 


CO 

a> 

Q. 
< / ) 
a> 

U) 

CO 

</) 

D 

0> 

c 

c 

<o 

JZ 

o 


T3 

i 

T5 

C 

<0 

CO 


(gp)  i|g  Jad  suojoqd 


Chart  164.  Comparison:  Efficiency  at  capacity,  orthogonal,  M  =  4,  soft  and  hard  decisions,  preamplified 
noncoherent,  heterodyne  or  preamplified  coherent,  homodyne  coherent,  photon  counting,  quantum,  or 
Gaussian  quantum. 
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Chart  165.  Comparison:  Efficiency  at  capacity,  orthogonal,  M  =  16,  soft  and  hard  decisions,  preamplified 
noncoherent,  heterodyne  or  preamplified  coherent,  homodyne  coherent,  photon  counting,  quantum,  or 
Gaussian  quantum. 
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Chart  166.  Comparison:  Efficiency  at  capacity,  orthogonal,  M  =  64,  soft  and  hard  decisions,  preamplified 
noncoherent,  heterodyne  or  preamplified  coherent,  homodyne  coherent,  photon  counting,  quantum,  or 
Gaussian  quantum. 
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Chart  167.  Comparison:  Efficiency  at  capacity,  orthogonal,  M  =  256,  soft  and  hard  decisions,  preamplified 
noncoherent,  heterodyne  or  preamplified  coherent,  homodyne  coherent,  photon  counting,  quantum,  or 
Gaussian  quantum. 
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Chart  168.  Comparison:  Efficiency  at  capacity,  OOK  and  orthogonal — M  (Mp\)  =  2,  3,  4,  8,  16, 32,  64,  128,  256,  512,  1024, 
soft  decision,  preamplified  noncoherent. 
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Chart  169.  Comparison:  Efficiency  at  capacity,  OOK  and  orthogonal — M  (Mp\)  =  2,  3, 4,  8,  16,  32,  64,  128,  256,  512,  1024, 
soft  decision — heterodyne  or  preamplified  coherent,  or  classical  Gaussian. 


Chart  170.  Comparison:  Efficiency  at  capacity,  OOK  and  orthogonal — M  (Mp\)  =  2,  3,  4,  8,  16,  32,  64,  128,  256,  512,  1024, 
soft  decision — photon  counting,  or  Gaussian  quantum. 
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soft  decision — quantum  ML,  or  Gaussian  quantum. 
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Chart  172.  Comparison:  Efficiency  at  capacity,  MPSK,  M  =  2,  soft  decision — preamplified  DPSK,  heterodyne  or 
preamplified  coherent,  homodyne  coherent,  photon  counting,  quantum,  or  Gaussian  quantum. 
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Chart  173.  Comparison:  Efficiency  at  capacity,  MPSK,  M  =  4,  soft  decision — preamplified  DQPSK,  heterodyne  or 
preamplified  coherent,  Bondurant  photon  counting,  quantum  HD,  Gaussian  coded,  or  Gaussian  quantum. 
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Chart  176.  Quantum  eigenvalues  (v*):  Orthogonal. 
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